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Existence, Uniqueness and Lipschitz Dependence for 
Pat lak- Keller- Segel and Navier-Stokes in IR 2 with Measure-valued 

Initial Data 

Jacob Bedrossian* and Nader Masmoudi^ 



Abstract 



We establish a new local well-posedness result in the space of finite Borel measures for mild 
solutions of the parabolic-elliptic Patlak-Keller-Segel (PKS) model of chemotactic aggregation 
in two dimensions. Our result only requires that the initial measure satisfy the necessary as- 
£Lh . sumption max a;g R2 fi({x}) < 8ir. This work improves the small-data results of Biler [4] and the 

existence results of Senba and Suzuki [35]. Our work is based on that of Gallagher and Gal- 
lay 33J, who prove the uniqueness and log-Lipschitz continuity of the solution map for the 2D 
Navier-Stokes equations (NSE) with measure-valued initial vorticity. We refine their techniques 
and present an alternative version of their proof which yields existence, uniqueness and Lipschitz 
continuity of the solution maps of both PKS and NSE. Many steps are more difficult for PKS 
than for NSE, particularly on the level of the linear estimates related to the self-similar spreading 
solutions. 

in ' 1 Introduction 
in 



The primary focus of this work is establishing a large-data local well-posedness result in the space 
of finite Borel measures for the parabolic-elliptic Patlak-Keller-Segel model in two dimensions: 

u t + V • (uVc) = An, 

-Ac = u. ' 



This system is generally considered the fundamental mathematical model for the study of aggrega- 
tion by chemotaxis of certain microorganisms [54} 145] |43| [42] . From now on we will refer to (11. ip 
as Patlak-Keller-Segel (PKS). The first equation describes the motion of the microorganism as a 
random walk with drift up the gradient of the chemo-attractant c. The second equation describes 
the production and (instantaneous) diffusion of the chemo-attractant. PKS and related variants 
have received considerable mathematical attention over the years, for example, see the review [43] 
or some of the following representative works [211^11521111^1531121113115311111912!. 

An important and well-known property of (jl.ip in two dimensions is that it is L 1 -critical: if 
u(t,x) is a solution to (II. ip then for all A 6 (0, oo), so is 

1 / t x 
u\(t,x) = _«!__ 
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It has been known for some time that (jl.ip possesses a critical mass: if H^olli ^ 8tt then classical 
solutions exist for all time (see e.g. [581 13 [HI O [8] ) an d if ||wo||i > 8tt then all classical solutions 
with finite second moment blow up in finite time [44 \ [50l [13] and are known to concentrate at least 
8tt mass into a single point at blow-up [58] (see also [40|, I57j). Another important property of (jl.ip 
that plays a decisive role in our work is the existence (and uniqueness) of self-similar spreading 
solutions for all mass a £ (0,87r). These are known to be global attractors for the dynamics if 
the total mass is less than 8ir [13] and for the purposes of our analysis, should be thought of as 
analogous to the Oseen vortices of the Navier-Stokes equations. When studied in higher dimensions, 
(jl.ip is supercritical and the dynamics are quite different, see for example [H [261 121 EHl EE] . Variants 
of (II. ip involving nonlinear diffusion which are critical in higher dimensions have also been studied 
[4"?l [HI EQl EH [lOl H6] (see also the related [39] ) . The parabolic-parabolic version of (jl.ip has also 
been analyzed in various contexts (see e.g. [57l [HI EU [53j [6] ) . We should also mention that variants 
of (jl.ip have been studied in the context of astrophysics (referred to as Smoluchowski-Poisson) as 
a simplified model for the collapse of overdamped self-gravitating particles undergoing Brownian 
motion (see e.g. [22j [23 [59] ) . 

The goal of the present work is to prove the most general local well-posedness result known for 
(jl.ip . We work with so-called mild solutions, motivated by similar notions used in fluid mechanics 
(other authors have also used this notion for (jl.ip ). See below for the full definition and discussion 
(Definition [1]), but the main idea is that these solutions satisfy (jl.ip as the integral equation 

u(t) = e* A /x - f e ('-*) A V • (u{s)Vc{s))ds, 
Jo 

and satisfy the optimal hypercontractive estimate sup tg ( t) 11^(^)114/3 < 00 (the self-similar 
spreading solutions show that the rate cannot generally be better). We show that there exists a 
unique mild solution to (jl.ip given initial data which is a positive, finite Borel measure \x that 
satisfies max^g^ n{{x}) < 8ir. Moreover, we also show that the solution map is locally Lipschitz 
continuous with respect to the total variation norm of the initial data. This is the most general well- 
posedness result possible in this space without considering weaker notions of solutions which can be 
extended past blow up (such solutions do exist [31], see below for a discussion). The mild solutions 
we construct are smooth for t > at least for some short time, which is not possible if there already 
exists a concentration with critical mass. Biler proved [I] using a contraction mapping argument 
that if the initial measure has a small atomic part then one can construct a unique mild solution. 
Senba and Suzuki [58] construct weak solutions only under the assumption max x . gK 2 /x({x}) < 8ir, 
however, their solutions are not a priori mild solutions and it is far from clear that more general 
solutions would necessarily agree with the mild solution. Hence, our proof also yields existence, 
which was to our knowledge open. 

Much of our work and motivation is a result of the similarities (jl.ip shares with the Navier-Stokes 
equations in vorticity-transport form 

f u t + V ± V- Vw = Au, 

Existence of mild solutions to (jl.ip with measure-valued initial data was proved earlier in |27l [37] 
and similar to (jl.ip . it is relatively easy to prove well-posedness if the initial data has only very 
small atoms. However, in [33] . Gallagher and Gallay proved that given arbitrary initial vorticity 
in the space of finite Borel measures, there is a unique mild solution to (jl.2p and the solution map 
is log-Lipschitz continuous with respect to the total variation norm (see also [361 [34] for a proof of 
uniqueness of the Oseen vortex with point measure initial data). 
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The proof of Gallagher and Gallay [33] uses an accurate approximate solution and an intelligent 
decomposition of the error between the approximate solution and the true solution, which is shown 
to be very small in an appropriate sense. A Gronwall-type estimate is used to prove that if two 
solutions have the same initial data then they must differ from the approximate solution in the same 
way and hence are equal. However, the argument is not quite a contraction mapping. Consequently, 
it requires the a priori existence of well-behaved mild solutions and yields log-Lipschitz dependence 
on initial data, but not Lipschitz. Our argument follows the same general principles set forth in 
|33| . however, we use a different decomposition which allows stronger results. In particular, unlike 
|33| , our argument is a true contraction mapping, and this allows us to prove existence of solutions 
as well as the Lipschitz continuity of the solution maps of both (jl.ip and (]1.2h (see Theorem [3]) . 

As in [33] , the approximate solution is constructed by guessing that near a large atomic concen- 
tration in the initial data, for short time, the solution to (jl.ip or (II. 2\i should look like a self-similar 
spreading solution, and elsewhere can be approximated by a linear evolution. In order to close a 
contraction mapping argument, some knowledge about the linearization around the approximate 
solution is necessary. A 'brute force' linear analysis is likely intractable, however, it turns out that 
knowing good spectral properties of all the well-separated pieces of the approximate solution is 
sufficient to close the argument. In particular, we need good spectral properties of the linearization 
of (jl.ip or (jl.2p around the self-similar spreading solutions. For NSE, the nonlinearity vanishes for 
radially symmetric data, which is why the Oseen vorticies are simply the self-similar solutions to 
the linear heat equation. Moreover, the linearization around the Oseen vortices is relatively easy 
to analyze, as it is a sum of the Fokker-Planck operator and an operator which is skew-symmetric 
in an appropriate Hilbert space. Nothing analogous to these properties hold in the case of the 
Patlak-Keller-Segel system: the self-similar spreading solutions solve a genuinely nonlinear elliptic 
system and the spectral properties of the linearization are far from trivial to analyze. One of the 
main tools for dealing with the linearization is the spectral gap recently obtained by J. Campos 
and J. Dolbeault [17] . As the work of Campos and Dolbeault has yet to appear in print and since 
their result is so centrally important to ours, we have included an independent proof of a version 
of their result which suffices for our work in Appendix $0 (obtained after hearing of their result). 
This spectral gap needs to be adapted to the spaces we are working in, similar to what is done in 
Gallay and Wayne [36J for NSE (see Proposition [5] below). 

An additional technicality that appears here is the fact that the velocity field for PKS is not 
divergence free. This makes most of the results of Carlen and Loss [19J on the fundamental solutions 
of linear advection-diffusion equations inapplicable. 

Global measure- valued solutions of (11. ip in the sense of Poupaud's weak solutions [55], which 
make sense even if there are mass concentrations, have been constructed by Dolbeault and Schmeiser 
in [31] by taking sequences of regularized problems and extracting a measure for u(t) along with 
an appropriate 'defect measure' to make sense of the nonlinear term. It appears that the resulting 
solution depends on the chosen regularization, as the formal dynamics derived by Velazquez [62} 
[63] are different from those constructed by Dolbeault and Schmeiser. Whether or not measure- 
valued solutions can be uniquely selected by physically or biologically relevant criteria remains an 
interesting open question. See [64] for some work in this direction. 

Let us end this introduction by summarizing some of the main difficulties compared to the study 
of the NSE in [33]: 

(a) For the PKS, the vector field is not divergence-free, hence we cannot use the results of [IH] 
on the pointwise decay and localization for the fundamental solution of the linear advection- 
diffusion equation. 

(b) For NSE, the existence of a mild solution with strong a priori estimates was already known, a 
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fact which Gallagher and Gallay exploit multiple times in combination with the results of |19j . 
For the PKS, the existence of mild solutions with such general initial data was not known. 

(c) The self-similar profiles G a are not linear in a, as is the case for the Oseen vorticies. The 
critical mass 8ir will appear in many places in our analysis. 

(d) The linear operators we have to deal with are harder than those that arise in the study of NSE. 
For NSE, due to the divergence-free property, these linear operators are a skew-symmetric 
perturbation of a Fokker-Planck operator. 

In order to overcome difficulties (a) and (b), we had to find a better decomposition of the error 
terms between the solution and the approximate solution. This gives better control of the error in 
norms which permit us to close a contraction mapping argument, allowing also the deduction of 
Lipschitz continuity of the solution map with respect to the initial data. To our knowledge, this is 
a new result even for NSE. To overcome difficulty (d), the spectral gap of Campos and Dolbeault 
and known spectral gap properties of general Fokker-Planck equations are important tools in the 
linear analysis. 



1.1 Results 

The precise notion of weak solution we are using is that of a mild solution, which are motivated 
by similar notions in fluid mechanics and have also been used previously in the study of PKS (e.g. 

EOT)- 

Definition 1 (Mild Solution). Given fj, £ Ai + (M. 2 ), we define u(t) to be a mild solution to ([Lip 
with initial data /u on [0, T) if the following are satisfied: 

(i) u(t) ^* ii as t \ 0, 

(ii) u(t) £ XT where 

XT = C w ([0,T];M + (m 2 ))n \u(t) : sup t 1/4 ||«(t)|| 4 /3 < oo I , (1.3) 

I te(o,T) J 

(iii) u(t) satisfies the following Duhamel integral equation for all t G (0, T) 

u{t) = e tA /i - f e (t - s)A V ■ (u{s)Vc{s))ds, (1.4) 
Jo 

with — Ac(s) = u(s) in the sense 

c(t, x) = j log \x - y\ u(t, y)dy. (1.5) 

Here C w ([0, T]; A4 + (M?)) is the space of u{t) which take values in finite non-negative Borel 
measures continuously in time with respect to the weak* topology. Recall the estimate on the heat 
kernel, 

l|e tA V/|| p <t-V2+VP-V9||/|| 9) (1. 6) 

which is a consequence of Young's inequality for convolutions. This estimate implies that (ii) ensures 
the Duhamel integral in (|1.4p makes sense for any mild solution. The heat kernel also satisfies the 
following precise estimate [37]: for all \i G Ai(M. 2 ) and p £ (1, oo], 

limsupt p\\e (J*\\ P < \\n\\pp, (1.7) 
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where || • || pp denotes the semi- norm on M(M 2 ) which measures the total variation of the atomic 
part: 

IImIIrp : = \rii x })\' 

{xm 2 :H{x})\>0} 

Estimate (jl.7p and related estimates play a key role in our analysis and the work of Gallagher and 
Gallay, as they show that size conditions for short time results should only depend on the atomic 
part of the initial data. Additional important facts about mild solutions are summarized in the 
following theorem. 

Theorem 1. (i) Let u(t) be any mild solution to PKS which exists on some time interval [0, T], 
T < oo. Then sup te (o,T) < oo for all p E [1, oo]. 

(ii) (Biler J^]/) There exists some eo > such that if Ho E A4+(M. 2 ) and satisfies 

limsupt 1/4 ||e* A // ||4/3 < e o, 
t\0 

then there exists a unique local-in-time mild solution to (jl.ip with initial data hq. 

Part (i), (to our knowledge new), shows that the condition sup te ? ,r) * 1/ ' 4 ]|ii(i) H4/3 < 00 is 
equivalent to the L°° hypercontractivity estimate u(t) ^t^ 1 (the proof shows that all such estimates 
are equivalent). Accordingly, standard parabolic theory implies that all mild solutions are smooth 
and strictly positive after t > until (potentially) critical mass concentration. Part (ii) is due to 
Biler [3] and combined with (jl.7p shows that given a measure with a sufficiently small atomic part, 
one can construct a unique mild solution local in time. Part (i) will play a role in the proof of the 
main results of the paper, although (ii) will not. 

We now state our main results. For PKS we prove the following existence and uniqueness 
theorem: 

Theorem 2. Let \i E .M+(R 2 ) with max xgK 2 fi({x}) < 8ir. Then there exists a unique, local-in-time 
mild solution u(t) to (jl.ip with initial data \x. 

As discussed above, our approach also yields the Lipschitz continuity of the solution maps 
for (II. ip and (|1.2p . Even for NSE, this is an improvement of the existing result of log-Lipschitz 
continuity, due to Gallagher and Gallay |33j. 

Theorem 3. The solution maps of both NSE and PKS in two dimensions are locally Lipschitz 
continuous with respect to the total variation norm. That is, for all fJ. 1 ,^ 2 E A4(M. 2 ) (in the case of 
PKS, we assume additionally jJ 1 £ A4 + (M. 2 ) and max//({x}) < Sir) with associated mild solutions 
w l {t),w 2 {t), there exists some constant Cl = Cl([i 1 ) > and T = T([i l ) > such that for all 
5 > sufficiently small, if 

II pt 1 - i^Wm < s, 

then 

sup (\\w\t) - w 2 (t)\\ L i + t^Ww 1 ® - w 2 (t)\U /3 ) < C L 5. 
fe(o,T) v ' 

Here, || • ||x denotes the total variation norm on finite Borel measures. 

Let us briefly discuss the energy structure of (jl.ip . which is important for characterizing the 
self-similar spreading solutions and for analyzing the global behavior, the former being crucially 
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important for our work. Formally, the Patlak-Keller-Segel model (jl.ip is a gradient flow in the I? 
Wasserstein metric for the free energy (see [8j [9] ) , 

F(u) = j u(x) log u(x)dx + — J J u(x)u(y) log \x — y\ dxdy. (1-8) 

In particular, if the initial data has finite free energy, then for reasonable notions of weak solution 
we have the energy dissipation inequality, 

T(u(t))+ [ u{s) |Vlogu(s) - Vc(s)| 2 dxds < -F(u(0)), 
J o 

for all t > until blow-up time. Using the sharp logarithmic Hardy-Littlewood-Sobolev inequality 
(see e.g. [18J) this implies global existence of any weak solution which has finite initial free energy 
provided that the total mass is less than 8ir [301 13] • The energy dissipation inequality is actually 
stronger in similarity variables: 

e = ^=, r = logt (1.9) 
and w(t, £) := tu(x,t). In these variables, (jl.ip becomes the following, 

w T + V • (wVc) = Aw + ±V • (£«;) ^ _ 

—Ac = w, 

which is formally a gradient flow for the self-similar free energy 

Q(w) = J «,(£) log w(£)d£ + w(0 |e| 2 d(L + ^j f w(Z)w(() log |£ - C| d£d(. (1.11) 

As the second moment is now part of the energy, uniform control on the entropy J w(£) \ogw(£)d£ 
from below can be obtained and the sharp logarithmic Hardy-Littlewood-Sobolev can then be used 
to show that any solution to (jl.lOp with finite self-similar free energy and mass strictly less than 8ir 
is uniformly bounded in time. In physical variables this is the optimal decay estimate u(t,x) < i -1 
as t — > oo. 

The free energy Q is important to characterize the self-similar solutions of (jl.ip . Biler et. al. 
show in [7J that for all a £ (0, 8tt) there exists a unique, radially symmetric self-similar solution with 
mass a, denoted here in self-similar variables by G a (£) (existence had been previously established in 
[5j [52] ) . These solutions will play the role that the Oseen vortices play in [33J as the approximation 
for the solution near the large atomic pieces of the initial data. The following proposition collects 
the important properties of the self-similar solutions, which show that in many ways they are 
qualitatively similar to the Gaussian Oseen vortices of the NSE. While these results are trivial for 
NSE, they are more difficult for PKS, due to the fully nonlinear nature of the self-similar solutions 
G a . Parts (i-iii) are not new, but we sketch some aspects of the proof in Appendix §|B] for the 
readers' convenience, as they are not all located in one place in the literature. Parts (iv) and (v) 
seem to be new and are both crucial in deducing the Lipschitz dependence in Theorem [3j Part (iv) 
is also necessary to show that many constants and linear estimates are uniform as a \ 0, which is 
necessary to prove Theorem [21 We should point out that the result of (v) depends on the spectral 
gap of J. Campos and J. Dolbeault. (see Appendix $C]for an independent proof). For any / £ L , 
/* denotes the Riesz symmetric decreasing re-arrangement (see [48] for more information on this 
symmetrization technique). In what follows we denote the polynomial weighted L 2 space, 

V(m) := {/ € IP : (0 m f(0 G IP} , 
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with the convention (£) := (1 + |£| 2 ) 1/2 . We also define Lg(m) = {/ G L p (m) : / /a!x = 0}. Note 
that for m > 2, L 2 (m) M> L 1 . For any m > 2 we have the following, which will be useful later 

t V4||i /(logti Lz£)|| 4/3 = n/(iogt,.)|| 4/3 < ||/(iogt,0IU 2M . (1-12) 

Now we may state Proposition [TJ 

Proposition 1 (Properties of the self-similar solutions). Let a G (0, 8tt). 

(i) There exists a stationary solution to (jl.lOp G a G C°° which is strictly positive, satisfies 
G a = G* a , ||G a ||i = ol and satisfies 

Ga{0 <a e" l?|2/4 for £ sufficiently large. 

(ii) G a (£) is the unique stationary solution of (jl.lOp with finite self-similar energy (jl.lip . More- 
over, G a is the unique minimizer of the self-similar free energy. 

(Hi) In physical variables, jG a * s the unique mild solution with finite self- similar free energy 

with initial data aS, where 5 denotes the Dirac delta mass. 

(iv) For a sufficiently small, the estimate \\G a \\Lv + HG^H^^) < Pim a holds for all 1 < p < oo 
and m > 2. 

(v) For all K < 8tt and for a, (3 < K, the estimate \\G a — G^\\lp + \\G a — Gp^jjtr^ < P ,m,K \ & — P\ 
holds for all 1 < p < oo and m > 4. 

Remark 1. Proposition Q] (iv) proves something slightly more precise: that \\G a — aGW^r^ + 
\\G a — aG\\ p < a 2 for a small, where = (47r) _1 / 2 e~^l / 4 is the standard Gaussian. 

Due to the a priori estimates and the uniqueness of G a , a compactness argument shows that if 
J u(t)dx = a, then lim^oo \\ju(logt, ^) — jG a (^)\\ p = Oforp G [l,oo]. These results are naturally 
analogous to the well-known results for the heat equation and for the 2D Navier-Stokes equations 
[36] . The spectral gap deduced by J. Campos and J. Dolbeault and the extension discussed below 
in Proposition [5] can also be used to deduce an exponential estimate on the rate of convergence, as 
in for example [361 H2] • 

Remark 2. An obvious question arises about whether or not Theorems [2] and [3] can be extended 
to more general models than PKS and NSE. If the nonlocal velocity law is a linear combination 
of the Biot-Savart law for NSE and the chemotactic gradient law for PKS then this generalization 
should be more or less straightforward since the G a will still be the self-similar solutions. However, 
if the velocity law is no longer homogeneous, for example if the equation for the chemo-attractant is 
replaced by — Ac + c = u or —V • (a(x)Vc) = u, then there are no longer exact self-similar solutions. 
If G a are still good short-time approximations for the evolution of atomic initial data, as should be 
the case in the examples just mentioned, then the stated results of Theorems [2] and [3] can likely be 
proved with similar arguments after some additional approximations. Such cases should also include 
models in which the chemo-attractant and/or the density u(t,x) is subjected to an external drift, 
provided the drift is sufficiently regular. If G a no longer provide a good short-time approximation 
to atomic initial data, more substantial changes would have to be made. 
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Notation and Conventions 



We denote the L p {dx) norms by \\u\\ p := ||ti||_Lp. If a measure other than Lebesgue measure is used 
to define the norm, this is denoted by L v {dp.) (note that this is different than the definition of the 
polynomial weighted space L 2 (m)). 

To avoid clutter in computations, function arguments (time and space) will be omitted whenever 
they are obvious from context. In formulas we will sometimes use the notation C(p,k, M, ..) to 
denote a generic constant, which may be different from line to line or even term to term in the 
same computation. Moreover, to further reduce clutter in formulas, we make very frequent use of 
the notation / < P) fc r .. g to denote / < C(p, k, ..)g. We will generally suppress the dependencies 
which are not relevant for the estimate at hand and simply write / < g. In most cases, universal 
constants from functional inequalities and parameters which are not important for the discussion 
are omitted. We will also usually suppress dependence from uniform estimates which have already 
been established, although we often alert the reader to the estimate being used. 

Self-similar solutions of mass a are denoted G a , and when the mass is given by oti, we will often 
shorten this to G{ . Similarly, the velocity field associated with the self similar solutions are denoted 
v Ga or v Gi . 

2 Preliminaries 

The following proposition collects the basic properties of the nonlocal velocity law of which 
are essentially analogous to the properties of the Biot-Savart law for NSE. 

Proposition 2 (Properties of the nonlocal velocity law). Define 

x 



Z7T \X\ 



Then 



(i) Let | = | — \ for some p G (1,2). Then, 

\\B*u\\ q <\\u\\ p . (2- 1 ) 

(ii) Let p G (1, oo). Then, 

IIV-B * u\\ p < \\u\\ p . 

Moreover, V • B * u = —u. 
(in) If u S L 2 (m) for some m £ (0, 1) or u G L^im) for some m G (1, 2). Then for all q G (2, oo) 

\\(0 m ~ 2 «B*u\\ q <\\u\\ LHm) . (2.2) 

2.1 Outline for the proof of Theorem [2] 

For e > chosen small later we define the decomposition of the initial data 

N 

fJ- = atj5 Zi + no, 

i=l 

for 5 Zi := 5(z — Zi), Zi G M 2 and a, > chosen such that ||/^o||pp < £■ If the measure [i contains only 
finitely many point masses then N is finite and independent of e for e sufficiently small. However, in 
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general there may be infinitely many point masses and in this case it is important to note that N is 
fixed large when e > is fixed small. Define the minimal distance between any two concentrations 
(which is generally forced small when e is chosen small): 

d : = min(|zj — Zj\) > 0. 

The goal of this decomposition is to construct an accurate approximate solution and use a pertur- 
bation argument to build a true solution which is very close to the approximate one. Analogously 
to [33], we construct a mild solution u(t) via a decomposition of the form 

u(t,x) = w (t,x) + Y^on-Wi \\ogt, -j= - j + -G a% y - j , (2.3) 

with the terms wo, Wi defined below. However, our definition of Wi is different than in |33j . 

In what follows we will explain the decomposition formally, assuming that we have a well- 
behaved mild solution already. In reality, we will construct this solution using the decomposition. 
The term wo is defined as the solution associated with the (approximately) non-atomic portion of 
the initial data, which formally satisfies the initial value problem 

d t w + V • (w v) = Aw , . 

wo(0) = Mo, 1 ' 

where still v = B * u is given by the nonlocal velocity law associated with the full solution. Since fio 
has a small atomic part, (II. 7ft suggests that for short time t 1//4 ||it;o(i)||4/3 will be small. Of course, 
it will take some work (Proposition [6]) to make this convincing, as v(t, x) is very singular at time 
zero. On the other hand, the part of the solution associated with the large atomic parts of the 
initial data is not small in any relevant sense, so further decomposition is necessary. Consider the 
solutions W{(t, x) of the advection-diffusion equation in physical variables: 

d t Wi + V • (wiv) = Awi 
Wi(0) = cti5 z 



In [33], the authors consider the difference between wi and the self-similar solution of mass ct:j 
centered at z%. This quantity turns out to be small as t~ 1 G ai ((x — Zi)t~ 1 / 2 ) is an accurate approx- 
imation for Wi for short time, however, Wi(t,x) — t~ 1 G cti ((x — Zi)t~ l l 2 ) proves difficult to correctly 
control in a contraction argument. Hence, we choose a different decomposition which still satisfies 

and while each Wi will be localized around Zi, Wi ^ W{ — G ai (although the proof will show they are 
close to being equal). In particular, (|2,5p cannot be decoupled into separate expressions for 
Applying (|2.5p to v = B * u implies a more precise PDE for wq: 



d tW0 + V-(w Ef=i ^o) + V • (w Ef=i aj^Vj (bgt, ^f)) = Aw 

wq(0) = Mo, 



(2.6) 



where % = B * uio and Vj = B * Wj . For future convenience define 

Vj(t,x) := OLj—pVj ( logt 
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We now turn to the definition of Wj for j > 1, which is somewhat more technical. For notational 
clarity define 



N 



(t,0 ■= ajVj(T,€- (zj - Zi)e r/2 ), 



2V 



(2.7) 



the velocity fields induced by the perturbations j ^ i in the coordinate system of Wi and the velocity 
fields induced by the self-similar solutions of the j 7^ i concentrations written in the coordinate 
system of Wi. Let 4>(x) be a smooth, non- negative, radially symmetric, non- increasing function such 
that 4>{x) = 1 for |x| < d/2 and cj>(x) = for |x| > 3d/ 4. We will define Wi to be a solution of the 
following: 

d T Wi + V • {wiV G ') + V • (Gm) 

+V • [wi Z^i (1 - <P(^ 2 + Zi- Zj ))v G * (£ - {z 3 - Zi )e-^ 2 ) 
+V • (Ei#j ftfiK ~ (*i ~ z^e-^me^v^) 

+V • (±G iV *) + V • (±GiV W <) + V • (±G^/ 2 v (e T , £e r/2 + *)) ( 2 ' 8 ) 

+V • (uJtJi + u^) + V • ( Wl e T / 2 v {e T ,£e T '/ 2 + Zi )) 
= Awi + |V • (£wi), 
lim r _ > _oo Wi(r) = 0. 

It is in the second and third line where our definition differs from [33] . Our definition more naturally 
treats the dangerous terms when making estimates, but destroys the advection-diffusion structure 
of (|2.8p and the coupling makes it trickier to prove that all mild solutions can be decomposed in this 
manner (proved below in Proposition [3]). We re-write the equations for the perturbations wo,Wi as 
the corresponding Duhamel integral equations. Define Sjv(i, s) to be the linear propagator for the 
PDE 



v Gj 



A/. 



Hence we may re-write (|2.6p as the formally equivalent Duhamel integral equation 



(2.9) 



wo(t) = £V(*,0)^o 



SN(t,s) 



N 



V- (u)o(s)£ (s)) + V- {w {s)J2 v j( s )) 

3=1 



ds. 



(2.10) 



The important properties of Sjv(i, s) are collected in Proposition [6] below. We now turn to the 
perturbations W{. Define the following linear operator, which is the linearization of the transport 
term around the self-similar solution, 

A a f :=V-(G a v) + V-(fv G «), 
v = B*f, 

and define the Fokker-Planck operator 



Lf := A/ + -V •(£/)■ 



(2.11) 
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Denote by T a (r) := e T ( L Aq ) the linear propagator for the PDE 

d T f = Lf- A a f. 

The important properties of T a are collected in Proposition [5] below. We may now write (|2.8p as 
the formally equivalent Duhamel integral equation 

Mr) = ~ r T at (r - r')(V • (wi £(1 - m T ' /2 + Zi- Zi))v G ^i - (z 3 - Zi )e~ T ' l 2 )))dr' 

(—GiV 9l ))d,T' - f T TaAr-r'XV-i-GiV^dr' 

^G ie r'/2^ e r'^ e r'/2 + z . )))dT f 
tti 

(aiWidi + WiV Wt ))dT 

{wie T ' l 2 v {e T ' ,ie T ' I 2 + Zi )))dr' . (2.12) 

The primary effort of proving Theorem [2] goes into showing that the system of integral equation 
(I2.10p . (|2.12p has a unique solution in the relevant spaces, which is done using a contraction mapping 
argument. The perturbations {wi}^ =0 are normed with M[w](t) defined as follows, which differs 
from the norm used in |33j by the presence of the constant Kq to be chosen later. Let 

M [w ](t) = sup s 1/4 ||fZ;o(s)||4/3> 

0<s<t 

and for 1 < i < N, 

Mi[Wi](t) = SUp \\Wi{ T )\\L\m), 
— oo<r<log(t) 

then define 

M[w](t) = max(K M [w](t),Mi[wi](t)) 

for some large constant Kq > 1 to be chosen later. We use Kq to enforce more control over wq than 
the other perturbations, which is very important for dealing with the potentially disruptive effect 
of wo on Wi, i > 1. 

By construction, the unique solution to the system (I2.10p . (l2.12|) can be re-constituted into a 
mild solution u(t) of (jl.ip via (|2.3p . However, it is not a priori clear that every mild solution can be 
represented as a solution to the system. This nontrivial fact is stated in the following proposition. 
The proof mainly depends on a compactness argument and that G a are the unique self-similar 
solutions, analogous to Proposition 4.5 in [33]. However, unlike [33], an additional step is required 
to construct Wj which satisfy (|2.12p since the cross-terms in (|2.12p couple all the Wj in a more subtle 
manner than in [33J. 

Proposition 3 (Equivalence of formulations). Suppose u(t) is a mild solution of (jl.ip . Then u(t) 
can be decomposed as in (I2.3P with wo,Wi satisfying the integral equations (|2.10p . (|2.12p . 

By Proposition [3j any mild solution must correspond to the unique solution of the system 
(I2.10p , (|2,12p , which would complete the proof of Theorem [2) 



- [ T T ai (r-r')(V- 

J — oo 

- [ T r ai (r-r')(V- 

J — oo 




T ai (T-T')(V- 
r ai (r-rO(V- 



11 



Remark 3. Alternatively, in order to complete the proof of Theorem [2] we could fall back to 
a proof which more closely matches Gallagher and Gallay and use their decomposition to show 
that any second mild solution must agree with the one constructed with the integral equations 
(|2.10p . (|2.12p . This should work, however, we prefer to give a more self-contained proof by going 
through Proposition El 

The rest of the paper is organized as follows. In Section §2.21 we state the main linear estimates 
which are required for the proof of both Theorem [2] and [3l In Section $3] we prove Theorem [2] 
in several steps. In Section $3] we establish Theorem [3j the proof of which is closely related to 
the main steps of Theorem [2j In Appendix $A] we establish the linear estimates stated in §2.21 
and in Appendix ^B] we sketch the proof of Proposition [TJ Finally in Appendix $0 we include an 
independent proof of a version of the spectral gap estimate due to Campos and Dolbeault. 

2.2 Requisite Linear Estimates 

We briefly recall some known properties of the linear propagator of the Fokker-Planck equation 
S(t) := e rL in L 2 (m), studied in [35J. The following proposition can also be found in [33] . 

Proposition 4 (Properties of S(r)). Fix m > 1. Then, 

(i) S(t) defines a strongly continuous semigroup on L 2 (m) and for all w G L 2 (m), 

\\S(T)w\\ L 2 {m) < \\w\\ L 2 {m) , \\VS(T)w\\ L 2 {m) < a ^y/ 2 IMU 2 (m)> ( 2 - 13 ) 

for all t > and where a(r) = 1 — e~ T . 

(ii) If m > 2 and w G Lg(m), then 

\\S(r)w\\ L 2 {m) <e- T l 2 \\w\\ L 2 {m) , Vr>0. (2.14) 

(Hi) If q G [1, 2] then for all w G L q (m) and r > 0, 

\\S(i~)w\\ L 2 {m) < -riTlklU«(m) ( 2 - 15 ) 

a(r)<? 2 

\\^S(T)w\\ L 2 {m) < rllHU«(m)- ( 2 - 16 ) 

a(r) i 

Note that 

V5(r) = e T/2 5(r)V. (2.17) 

The following proposition is of crucial importance. It is the analogue to Proposition 4.6 in 
[33], but the proof deviates in several key places due to the different nature of the linear operator. 
Indeed, Recall that analyzing the spectral properties of the linearization around G a is more difficult 
for PKS than for NSE, as the operator A a is not skew-symmetric in any relevant Hilbert space. As 
mentioned previously, the key tool used is the spectral gap recently obtained by J. Campos and 
J. Dolbeault [IT] . This spectral gap must be adapted to the polynomial weighted spaces L 2 (m), a 
procedure analogous to what is done in [33], which we carry out in Appendix §Al 

Proposition 5. Fix a G (0, 8n) and m > 2. 
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(i) Ta( T ) defines a strongly continuous semigroup which is bounded on L (m) and satisfies 

\\r a (r)f\\ L 2 im) < a \\f\\ L 2 (m) , (2.18) 

l|Vr Q (r)/|| i2(m) < Q a(T)- l / 2 \\f\\ L 2 {m) , (2.19) 

where a(r) := 1 — e~ T . 

(ii) For some v = v{a) G (0, 1/2) which depends on a and for all f G Lq(tti), 

\\T a (r)f\\ L 2 (m) < a e- VT \\f\\ L 2 (m y (2.20) 

(Hi) If q G (1,2] then 7^(r)V is a bounded operator from L q (m) to L^{m) and there exists exists a 
v G (0, 1/2) (the same v as in (ii)) such that, 

e~ VT 

\\r a (r)Vf\\ LHm) < a -—-ll/ll^j. (2.21) 

Though v and all of the implicit constants depend on a, as a \ 0, v ~ 1/2 and the constants 
are uniformly bounded by Proposition^ (iv) as T a {j) can be treated as a perturbation of S(t) (see 
Remark^ in Appendix fC]). 

The following proposition is the analogue of Proposition 4.3 in [33], but the proof must deviate 
from the corresponding one for NSE in a non-trivial manner, as the underlying linear operator 
no longer has as nice of a structure. Our proof uses spectral properties of linear Fokker-Planck 
equations with general confining potentials. The proof is carried out in Appendix 



Proposition 6. There exists some to sufficiently small (depending on d and M pp ) such that 

(i) SN(t,s) defines a strongly continuous semigroup on A4(M?) and for all p G [l,oo] and v 6 
A4(R 2 ) we have 

\\S N (t,s) 

v I \lp ^ Tj— u_i/ p IMIa4> 0<s<t<s + to. (2.22) 



(ii) For allp £ (l,oo] and v 6 M + (M 2 ), 



limsupt 1 - 1 /!^,^!!,, < |M|pp. (2.23) 

t\o 

(Hi) There exists some Ao £ (0,1/2) independent of e such that the following holds: for allp G 
[l,oo], for all 7 G (0, A ) and for all w G L 1 (R 2 ), 



I / t \ 7+1/2- Ao 

(t - sf/ 2 ' l /P \s, 



\S N (t,s)Vw\\ p < Pi7 _, a/9 _ 1/w ( - ) Hli> 0<s<t<s + t . (2.24) 



All of the implicit constants above are independent ofe,N and d. Moreover, it will suffice to choose 
to such that 

t < d 2 min(l,iT), 
for some K which is independent of e, N and d. 
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3 Existence and Uniqueness: Proof of Theorem [2] 



We proceed in several steps. First we prove the contraction mapping argument which establishes 
the existence and uniqueness of solutions to the integral equations (]2.10p and (I2.12P which is the 
core of the proof. Next we establish Theorem Q] (i) which is necessary to establish Proposition [3l 
which is carried out last. Finally we briefly summarize the full argument at the end of the section. 

3.1 Contraction Mapping 

We will construct our solution to (|2.10p and (|2.12p in the following ball (for e > and T > to be 
chosen small later): define p := {pi} i=0 , 



Note that the ball is centered around Sjv(i, 0)po, although given Proposition [6] (ii), this is a minor 
detail. For any {u>i}^L , the corresponding u(t,x) constructed by (|2.3|) will be in xt (but is not 
small due to the presence of the large atomic pieces). It might be useful to bear in mind that the 
approximate solution we are perturbing around is, 



although we will not make explicit note of this in the remainder of the paper. 

Let w = F[p] be the nonlinear solution map which takes p to w := {u>i}^Lo defined by the 
following procedure. In what follows, define 



B e>T = {p(t) : M[p - S N {t,0)p ,p l > 1 ](T) < e} . 




B * p , 
B* Pj , 





Given p, we define wq by 



(t) = S N (t, 0)p - [ S N (t 



A' 



Jo 



s) V-(p (s)v p (s)) + V-(p (s)Y,^vf(s)) ds. 



i=l 
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Similarly we define u>i by 

Mr) = ~ 1^ T ai (r - r')(V • ( Pl ^(1 - <f>(^'/ 2 + z t - Zj ))v G i (£ - (zj - z % )e~ r ' / 2 )))dr' 
J -°° i#< 

- [ T T ai {r-T')(V -{-G^dr 1 - f r ai (r-r')(V-(-G^))dr' 

J — oo J —oo 

- ^ r Ql (T-T / )(V-(-G^'/2^ (e -' )ee -72 +Zi)))dr / 
J — oo &i 

T ai (r - r')(V • (^a^f + Pi v Ri ))dT' 
T ai (r - r')(V • ( Pl e T '/ 2 v p (e T ',Ze T/2 + Zi)))dr' . 



By definition, a fixed point of F corresponds to a solution of the system (|2.10p . (|2.12p . which in turn 
corresponds to a mild solution of 

The application of the contraction mapping theorem requires that F maps B ty T to itself and 
that F defines a locally Lipschitz mapping, which we prove in separate propositions. Note that 
linear terms in the definition of F can be treated essentially the same in the two propositions, and 
it is in the linear terms that our arguments significantly differ from |33| . 

Proposition 7. For e > sufficiently small, T > can be chosen sufficiently small such that 
F : B e>T -> B e , T - 

Proof. First we control wq to show that for T chosen sufficiently small (depending on e) 

K M [w - S N (t,0)p ](t) < if!if Mo[p](t)M[p](i) < K x (M[p)(t)) 2 , (3.1) 

for some constant K\ > which is independent of e. This is essentially the analogue of Proposition 
5.1 in [33 J, and we approach it in a similar way. Indeed, 



Kot^WMt) - S N (t,0)p \\ 4/3 < Kot 1 ^ / \\S N (t,s)V ■ (po(s)v p (s))\U /3 ds 



K t 1/4 \\S N (t,s)V-(p (s)J2<Xitf(s))\U /3 ds. (3.2 



To control the first term: by Proposition [6l Holder's inequality and (|2.ip we have for t sufficiently 
small (so that Proposition [6] holds), 



\\S N {t,s)V ■ (po(sH(s))\\ i/3 < 1 (^) 7+V2 A ° UpoWSgOOIIi 



1 


(t 


- S) 3 / 4 






(< 






1 



7+1/2-Ao 
7+1/2-Ao 
7+1/2-Aq 



^ 1 - 1 \\P0(S)\U/3\K(S)\U 
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Hence, 

i^ 1/4 f \\S N (t,s)V ■ (p (sH(s))\U/sds < Kot^j* 1 (;) 7+1/2 A ° \\Po(s)\\l /3 ds 

< Kof* (M M{ t)f ^ ^-1^ (J) 7 " 72 " -Lds 

<i^ 1/4 (M [p ](t)) 2 . 
To control the second term in (|3.2[) we proceed similarly (again for t sufficiently small), 

Kot 1 / 4 / \\S N (t,s)V ■ (p (s)J2^j^))h/3ds < 
J o j=1 

N ft 1 / + \7+l/2-Ao 



^2ajK t 1/4 y~ J IIPo(s)||4/3l|Pi(s)IU/3^- 



Recalling (fTTT^l . 

i^ 1/4 / ||5jv(t,«)V-(po( 8 )X;«iSfW)IU/3d«< 
J o j=i 

N rt I / + \ 7+1/ 2 - Ao 



Putting the estimates together proves (|3.ip . 

The significantly more delicate challenge is controlling Wi for 1 < j < N, which is stated in the 
following lemma. 

Lemma 1. For all e > the following estimate holds for i £ {1, N}: 

Mi[w){t) < 62(f) + v(t)M[p](t) + K 2 M [p ](t) + K 3 Mi[p](t)M[p](t), 
< 6 2 (t)+ V (t)M[p](t) + ^M[p ](t) +K 3 M i [p](t)M[p}(t), 

where 62(f) and rj(t) depend on e but go to zero as t \ while K 2 and K3 are independent of e. 
Proof. Write (|2J2j) as 



l (r) = Y J F l , k (r) 



k=l 
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where 



p iM = E / T T ^( T - r ') V • (- G ^ Gj (^ " & - Zi)e- T ' /2 ))dr' 
F h2 {r) = J2 r T ^ T ~ r ')V • {-G^v^i - (zj - Zl )e- T 'l 2 ))dT' 

... J — oo a i 

Fi, 3 (r)= f T ai (r -r')V ■ (-G ie T '/ 2 v p (e T ' ,^e T '/ 2 + z t ))dr' 

J — oo 

FiAr) = F$( T ) + if 4 } (r) = f T ai (r - r')V • ( Pi £(1 - <f>(&'/ 2 + z t - Zj ))v^(C - { Zj - z^' l 2 ))dr' 

J -°° j& 

+ [ T T ai (r-T')V-(Y,^P^ + ^i-^ T 'hm T ' /2 )v Gi mdr' 

N 

j=l J-oc 

F ifi {r)= f T at (T-T')V-( Pl e T '/ 2 v p Q (e T ',^'/ 2 + z t ))dT'. 

J —oo 

The first is controlled analogously to the corresponding term in [33]. Using (I2.2ip . 

-u(t-t') r< 



l^i(r)|U 2M < Y, Q [ r 3)i/ 2 II^^ - & - Zi)e- T ' /2 )(0 m ^hdT> 
£ v ' J-oo a{r ~ t ') i/2 



with an implicit constant which is independent of e and N. The last inequality follows from Proposi- 
tion[T](iv) which implies {£,)~ 1 v Gj '(£ — (zj — Zj)e~ T//2 ) < C(aj)e T//2 , where the constants are uniformly 
summable in N. To control the next term we also proceed analogously to |33j . 

m,2\\m m) £5>/ T llff ft - (** - Zi^M^hdT*. 

j^i oo V 7 

Continuing with q G (2/V, oo) and 7 G (2/g, 1), 

n^(£ - (z, - , l )^ T ' /2 )(o m ^iiL 2 < iko 7 ~%ii 3 ii<£ - - z^'/^itr^w^- 

Cxi £*i 'J -2 



The first factor can be controlled via the weighted estimate on the nonlocal velocity law (|2.2p which 
implies (since m > 7): 

ll(0 T "%ll«<l|Pilk 2 ( 7 )<l|PilbM- 

The second factor is controlled by the localization of G a given in Proposition [U In particular, since 
d > we have, 

ii<* - (*,- - zo^> 2 Mfr-ii j* < d 
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Therefore, 



L 2 (m) 



-v{t-t') 



< 



M[ P W)e 



a(r - r') 1 / 2 



^^(7/2-1/9) 



Now we confront the next term using the LP estimate on the nonlocal velocity law (|2.ip , (|2.2ip and 
Proposition [Q 



ll-^,3||z,2(m) 



-v(t-t') 



< 



M [p](e 



a ( r _ T /)l/2"- -0 
e -^(r-r') 



dr' 



a(r - r') 1 / 2 



< ^ 2 M [p](e T ; 



Note that is independent of d and e, since by Proposition Q] (iv), for a — >• 0, ||(£) 4 %Hl4 remains 
bounded. Now we turn to ^4, which is an important difference between the work here and 
Dealing with the first term: 



ll^i,4 IU 2 (m) 



-u(t-t') 



o(r - r') 1 / 2 



||p, £(1 - 4>^ T ' 12 +Zi- Zj ))v G > (£ + {zi - z )e^' 2 ))\\ L , {m) dr'. 



Zi Zj 



> 



By the definition of the cut off <f>, the integrand of the L 2 norm is only non-zero if ^e T //2 + 
d/2, which of course is equivalent to £ + (z^ — Zj)e~ T> I 2 > e~ T '/ 2 d/2. Since v Gj decays like 
{ aj /2v)\t\-\ 

fT e -v{T-Tl) 



I 77^(1) II , < 



liT-T') 1 / 2 



\\ Pi £(1 - <P(^'/ 2 + Zi - Zj))v G '{Z + (Zi - z 3 )e^'' 2 ))\\ L , {m) dr' 



< 



3& 



a(r-T') 1 / 2 \ d ^ aj ) WP^ T ')\\^{m)dr' 

<^ 2 M[ Pl ](en. 
Now we consider the second term, which is concentrated around Zi, 

fT e -u{T-r') 



lKi,4 \\L 2 (m) 



^.no air - r')V 2 ^ ^ 



Since 1^(0 1 ) | is bounded uniformly (for a < 8tt obviously) by Proposition [IJ 

-\\ Pj (C - (zj - z t )e-^/ 2 me r/2 )v G ^0(O m h < \\PM ~ & ~ z^'/ 2 )^/ 2 )^ 1 h- 
Due to the cut-off, the integrand is only supported where |£| < 3e _T / 2 d/4 and this implies that 

i - ( Zj - Zi)e- T '/ 2 > e- T '' 2 d/A. 
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Hence, 



(0m-l ^ e -r'/2^m-l 



(£ - (Zj - z i )e-^'/2)m ~ (de-T'/ 2 Y 

< 1 



Therefore we can translate the coordinate system in the L 2 (m) norm and we get from the above: 

^ e -r'/2) ll ^ llL2 ( m )- 



< 



Putting the previous two estimates together we ultimately have 

||if 4 V) + F${T)\\ L * {m) < + j^j^j MH(e^), 

where the implicit constant does not depend on e or d. Estimating the first nonlinear term, for 
1< p < 2 using ([2T2TD 

._j ./-oo a \J T ) 

The norm can be estimated with Holder's inequality and the LP estimate for the nonlocal velocity 
law (|2.ip (and that L 2 (m) injects into LP ', p < 2 since m > 2): 

Oj - Zi)e~ T ' /2 ,T)pi(T')\\ L p {m) < ||(0 m pi||2||^ P ||2p/(2-p) < l|Pi|U 2 (m)INIp i$ IIPi||L 2 (m)ll/5illL2( m )- 

Hence, 

[|F f)5 || xa(m) <M i W(e T )M[p](e T ), 
with an implicit constant which is independent of e. The last nonlinear term we deal with similarly, 



< 



By definition, since t = e T , and using the L p estimate for the nonlocal velocity law (|2.ip we have 

5 P Wll4 



| e ^P (e -' ) ^'/2 + ^ )||4=t l/4 || ~ r 



<t 1 /4|| /0o || 4/3 = Mo[p]W . 

Therefore, 

\\F lfi (T)\\ LHm) <M [p}(e T )M[ Pl }(e T ), 
with an implicit constant which is independent of e. This completes the estimate of □ 
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In order to prove that F : B e j< — > B e j> we need to first fix e small, then T small and Kq large. 
The reader is advised to note that fixing e small in turn generally fixes N large and d small. Now, 
if we fix e, then restrict T such that 82(f) < e/4, rj(t) < e/4, and Kq such that Kq > AK2 then 

K M [w - S N (t,0)p ](t) < KiK M [p}M[p] < K t e 2 , 

and 

M[wi](t) < \ + \M[p\ + X -M[p\ + K 3 Mi\p]M\p] 
< 3e/4 + K 3 e 2 . 

Finally, the result follows by choosing e small. □ 

Proving that F is a contraction does not pose any significant new challenges. The linear terms, 
which were the most difficult to deal with in Proposition [7J are treated exactly the same. The only 
variation is in the treatment of nonlinear terms, but these do not pose a significant issue and can 
be dealt with as in [33]. Hence we only sketch the proof. 

Proposition 8. For e > sufficiently small, T can be chosen sufficiently small such that F is a 
contraction on B e T- That is, if p l ,p 2 E B e T and w 1 = F(p 1 ),w 2 = F(p 2 ) then, 

M[w x -w\t) ^^Mlp 1 - p 2 ](t). 
Proof. We first estimate Mq[w 1 — w 2 ](f). By definition, keeping notation analogous with above, 

wl(t)-w 2 (t) = - I S N (t,s)V ■ {i p 1 (s)p 1 (s)-v p 2 (s)p 2 (s)) ds 
Jo 

-J>, / S N (t,s)V ■{vf{s)pl{ S )-vf{s)pl(s))ds. (3.3) 
3=1 Jo 

All the terms are dealt with essentially the same, but consider the interactions of po with pi, i > 1. 
Write 

v^pl ~ vfpl = {vf - v?)pl + v?(pl - pi). (3.4) 
Estimating the first set of terms, using (|2.24p . Holder's inequality, (|2.1|) and (j 1 . 1 2 1) . 



N ft -, / / \7+l/2-Ao 

^E^S^ {t _ 8)3 /4 {-) \K(s)-V?(s)\\4pl(8)\U /3 ds 

ft 1 / f \ 7+1/2- A -, 

< m[„ - *HWVW"J (-) ^ 
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Similarly, 



N 



t 1 ^^ S N (t,s)V- ^pa^(s) (pi( S )-p2( S ))U|| 4/3 

N ft i /.x 7 +l/2-A 

^E* 1/4 «i/ o {t _ s)3/4 (-J ll«fWll4l|pSW-p8Wllv8* 

/■< i / + \ 7+1/2- A -, 

The terms involving Vq' are treated similarly (easier in fact) so we omit the details. Using these 
estimates together with (|3.3|) implies that there exists some constant K$ independent of e such that 
for all e > sufficiently small, we can choose T sufficiently small such that if t G (0,T), 

K M [wl - w 2 ](t) < K 5 M[ P1 - pa] (t) (AoMofp 1 ]^) + M[p 2 ](t)) 

< K 5 M[ Pl - p 2 ](t) (M[p l ){t) + M[p 2 ](t)) ., (3.5) 

Now we turn to the contraction estimate on the perturbations around the self-similar solutions. 
Similar to [33], define Gi : k = F^ k — Ff k , where F- k is defined as in Proposition [7] corresponding to 
p 3 . The source term satisfies Gn = 0, whereas the proof of Lemma Q] immediately implies that the 
linear terms satisfy 

G 2 , k + G 3:k + G 4 , fc < rj^Mlp 1 - p\t) + ^Mlp 1 - p 2 ](t), (3.6) 

where rj(t) and K 2 are the same as those in Lemma[TJ Finally, the nonlinear terms G^k, G$ t k can be 
treated easily by combining (|3.4p with the arguments of Lemma [T] to prove that there exists some 
Kq independent of e such that 

G 4 , fc + G 5 , fe < K 6 (M[p l ]{t) + M[p 2 ](t)) M[p l - p 2 ](t). (3.7) 

Together, (|3.5p . f|3.6[) and (|3.7p imply Proposition [8] by first choosing e sufficiently small: 

£< l^ mm (i<i)< 

choose Kq > 8K2 and then choose T such that (|2.24j) holds and that ry(i) < | (the parameters of 
course should also be chosen such that Proposition [7] holds). □ 



3.2 Proof of Theorem □] (%) 

We now prove Theorem [T] (%), which is an important property of mild solutions and also plays a 
necessary role in the proof of Proposition [3l 

Proof. By standard theory on the continuation of classical solutions, it suffices to assume T > is 
sufficiently small. We first prove that the ||it(£)|| 4 /3 < t -1 / 4 estimate can be bootstrapped up to 
the estimate ||w(£)||3 < i~ 2 / 3 . Then we show separately that this additional estimate implies the 
L°° estimate. This kind of two step bootstrap approach is common when applying similar methods 
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HH S3 13l [TOl [T31 ESI [3l [2] . It turns out to be a little easier to prove the L 3 estimate in self-similar 
variables. Hence define, 

w(r,£) =tu(t,x), r = logi, £ = ~/|- 

The hypercontractive estimates < i p are all equivalent to ||io(t)|| p < 1, which is one of the 

reasons the self-similar variables simplify the argument. We use an argument which takes advantage 
of the a priori control on the vertical distribution of mass implied by ||^(t) H4/3 < 1. Indeed, 
estimates on the vertical distribution of mass have long been known to be a key controlling quantity 
for PKS (see e.g. [1U [131 E]) so it is natural that such control also produces hypercontractive 
estimates, as already seen in [13j . Let k > 1 be some constant which will be chosen later and define 
w k = (w — fe)+. Let To G (—00, log(T) — 1) be arbitrary and define p(r) = 4/3 + (5/3) (t — To). Note 
that while p varies with t, it lies in p £ [4/3,3] and p = 5/3, so for making most estimates we can 
treat p as basically constant. We will show that there exists some constant C\ independent of to 
such that 



K1+to)||I = |K1+to)CS<C' 1 , (3.8) 



which proves the desired claim. For r E (to, to + 1) we now compute the following (defining 
—Ac = w), 

^- J w k {Tf^di = p J w p k \ogw k d^ + p J wl- l {Lw-V -{wVc))dt 

4(p - 1) 



p I w p k \ogw k d£ 



P 



Vw 



p/2 



d( 



+ 



J w p ^V ■ (£w)d£ -p j w p ~ l V ■ (wVc)d£. 



Using w l k w = w l k +1 + kw l k and —Ac = w, the above expands into 



^ J w k { T y^di=p J w i\ogw k dt-^—^ J | v < /2 |V + (p-i) J w p k +1 dti 



+ C(k,p) I w p d£ + C{k,p) j w p ~ l d^ 

for some constants that depend on k and p that we will not need the precise values of. Since 
p = 5/3 > and logitf^ < w k , 



d_ 
~d7 



w k (r) p ^dC < 



4(p ~ 1) 
p 



Vw 



p/2 



dC+(p+i) / w p+1 dH + C(k,p) I v%d£ + C(k,p) I w p 



We may interpolate all of the lower order terms between L p+1 and L 1 and use weighted Young's 
inequality to deduce 



1 



C(k,p)\\w k \\ p < C(k,p)\\w k \\ p ^ \\w k \\f < ^\\w k \\ p p X\ + C(k,p)\\w k \\i, 

where the constant in the last inequality is different than the first, but we do not need to track such 
details. A similar inequality holds for p — 1 and hence, defining M = 



d_ 



J w k {r) p ^di < 



4(p-l) 
p 



Vw 



p/2 



d£ + (p+l) I w p+1 d£ + C(k,p)M. 
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Similar to [15] . we apply the Gagliardo-Nirenberg inequality 

\\w k \f p Xl<K(p)\\ Wk \\ 



2 



to the first term, which implies, 



1 1 < - %£ "j^r + 6. + DIMS! + 

Now we use uniform vertical control imposed by the estimate 1 1 1 1 4/3 < 1, which implies, 

Applying this to the time evolution of ||u7fc||p, we have that for for some constant which is uniformly 
bounded for p on r G [to, to + 1], 

^ 1 u> fc (r)fWde < (p+ 1 - C^A: 1 / 3 ) |K||^ + C(k,p)M. 
Hence, for k chosen sufficiently large we have, 

A J Wk (rY*M < -\\w k \\ p p Xl + C(k,p)M. 
Since ||iOfc||p < H^fcllp+i + ll^fclli — ll^fcllp+i + M we finally have 



A J Wk { T y^dt<-\\w k \\i + c* 



where C* is some constant which depends only on M, asp € [4/3, 3] and k has been fixed. Integrating 
implies 

11^(^ + 1)111 < max (||tt>fc(TD)[|^3, C*J < max M|io(r ) ll^g, C*^ . 

Using now the inequality 1 1 1 1 3 < H^fclb + & 2 |M|i we finally get (|3.8|) . Notice that the a priori 
estimates on w(t) imply that ||-B*u>(r)||oo ^ 1, which in physical variables is equivalent to ||w(t)||oo ^ 
/ 1 2 . 

To bootstrap the L 3 estimate to L°° we return to the original variables, although the reader may 
wish to note the parallel between the following argument and the one just finished. Let t k = 2~ k 
and consider the dyadic intervals [t k ,t k —x\. In the following computations it is important to keep 
in mind that t k _\ = 2t k = 4£ k +i, and are hence all comparable. On each dyadic interval, 



u 



(t) = e {t - tk)A u(t k ) - / e {t ~ s)A V ■ (u(s)v(s))ds 



h: 
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By the L 4 / 3 and L 3 bounds on u, (fL6|) and |Ki)|U ^ t' 1 / 2 , 

Ht)\\oc <(t~ t fc )- 3/4 ||n(t fc )ll4/3 + f ||e (i ~ s)A V • (u(s)v(s))\\oods 

Jtk 



<q 1/4 (t-t k )-v 4 + / (t- S )- 5 / 6 ||n( S ) V ( S )|| 3 d S 
Jtk 

< t~ 1/4 (t - t k y 3 / 4 + fit - s r 5 / 6 |K s )|| 3 |K s )||oc^ 

Jtk 

< q l/ \t - t k )~^ + f\t - S )- 5 /6 s -2/3 s -l/2 ds 

<t- k i i\t-t k )-^ + ti 7 '\t-t k )^ 

Hence this implies ||u(tfc-i)||oo ~ tZ 1 ~ t k ~i- Re-doing the above computation using this informa- 
tion to deal with the first term (and the maximum principle ||e tA /|| 00 < ||/||oo) we see that indeed 

N^llooSt- 1 . " □ 

3.3 Proof of Proposition [3] 

In this section we prove the equivalence of the integral equations (|2.10p . (|2.12p with general mild 
solutions. Let u(t, x) be a mild solution of (jl.ip with initial data /i with associated nonlocal velocity 
v(t) = B * u(t). By Theorem [TJ u(t,x) necessarily satisfies the following a priori estimates for all 
p £ [1, oo] and q £ (2, oo], 

Nt) Hp < t 1/p ~\ \\v(t)\\ q = \\B * u(t)\\ q < iH. (3.9) 
Suppose for any e > we write 

N 

A* = HQ + / J a j S z i , 
3=1 

with as always ||/io||pp < e - Then we may define wq as the mild solution to (|2.4p and define Wi, 
1 < i < N as mild solutions to 

d t Wi + V • (wiv) = Awi 
Wi(0) = ai6 Zi . 

By definition, u(t, x) = uio-\-Yli=i w i- Note that since the W{ are all non-negative, each is guaranteed 
to still satisfy (|3.9|) (for NSE this is guaranteed by the results of Carlen and Loss [12] ; see also [33J ) . 
For alH £ {1, N}, we write Wi in self-similar coordinates (without re-naming), r = logt, £ = 
which then satisfy 

d T Wi + V • (viWi) + V • (RiWi) = Lwi, 
for Ri given as follows: define V{ = B * W{ then as in [33], 

N 

R l (r,O = e T/2 v (e T ^e T / 2 + z i )+ £ Vj {r,t - ( Zj - z^' 2 ). (3.10) 
The steps to proving Proposition [3] are the following: 



24 



(i) Show that ^(t, £) — > G a (£) in L 2 {m) for m > 2 as r — >• — oo. 

(ii) Use (i) to construct uij which satisfy (|2.12p . 

(iii) Show that wq satisfies (|2. 10|) . 

Once we have completed (i) and (ii), (iii) follows from the continuity properties of Sjsr(t, s), indeed, 
once we have constructed suitable Wi, we may write 



ttfo(i) = S(t,t )w (t Q ) - S N (t 

Jto 



N 

V • (wo(s)vo(s)) + V • (w (s)^2vj(s)) 



ds, 



and pass to the limit to \ using the Duhamel formula (IA.3j) in §A,2I and Proposition [H 

Now we concentrate on the more involved procedure of proving (i) and (ii). Part (i) uses an 
energy/compactness argument analogous to the approach of Gallagher and Gallay [33] , The idea is 
as follows: a compactness argument shows that iUj(r) is precompact in L 2 (m) as r — > — oo and the 
uniqueness properties of the self-similar solution stated in Proposition [T] will imply that the a-limit 
set can only consist of {G ai }. The first lemma is the compactness. 

Lemma 2. For all i £ {1, N}, {wi(r)} is precompact in L 2 {m) dsr-> — oo. 



Proof. The result of Carlen and Loss [19] for compressible advection-diffusion equations together 
with the a priori estimate on the velocity field given in (|3.9p imply the Gaussian localization estimate 



Wi (r,Oe^ d£<l, (3.11) 

for some 7 > 0. Combined with the a priori estimates (|3.9p . Wi{r) is then uniformly bounded in 
L 2 (m) for all m. As H 1 (m + 1) r )' > L 2 (m) by the Rellich-Khondrashov embedding theorem, it 
suffices to prove that Vwj(r) is uniformly bounded in L 2 (m) for all m > 2, for which we proceed 
similar to what is done to prove analogous statements in [33, 36J, with the necessary alterations to 
deal with the divergence of the velocity field. Write W{{t) in integral form: for some —00 < tq < 
log(T), 

Wi(r) = S(t - T )wi(T ) - S(t - s)V • (vi(s)wi(s) + Ri(s)wi(s))ds. 



Using (|2.13p . for some p £ (1, 2) we have, 



|V^(r)|| i2(m) < ^ T0 J h ^ + f T o(T _ s )-Vp|| V . ( Vl ( s ) Wi ( s ) + RitfwiisMvMda 
& \T i~o) J to 

< llw f To)h *$ + [ T a(r - S )- 1/p ||K( S ) + i2i(a)) • Vw^h^ds 
a(r - sr 1/p |Ms)(V • Ui(«) + V • Ri(s))\\ LP{m) ds. 

TO 



The second term can be controlled as in [33J: 

\\(vi(s) + Ri(s)) ■ Vwi(s)\\ LP ( m) < ||Vw i (s)|| £ ,2( m )||u i (s) + Ri(s)\\jp_. 



2-P 



The latter factor is bounded as follows, recalling the definition of Ri (|3.10p . For j G {1, .., , N} using 
(|2.1|) and Holder's inequality, 

IKOOII^il ^ lki(«)|| P <m \\Wj(s)\\ L 2 {m) < m 1, 

2-p 
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and for the approximately non-atomic part, using the a priori estimate (|3,9p 
\\e T / 2 v (e T ,^ 2 + zi)\\j^ = e T ^-^\\v (e T ,x)\\^ 

p— 2 p-2 

= t^-^\\v (t,x)\\^p_ < 1. 

p-2 

We now turn to the second term, which involves the divergence of the velocity fields. Recall that 
for PKS, the nonlocal law is given by v = V(— A) _1 u> and hence for all j G {1, ...,iV}, 

IK(t)V • Vj(r,£- (zj - Zi)e' T/2 )\\ LP{m) = \\wi(T)wj(T, £ - (zj - ^)e~ r/2 )|| L p (m) 

< IK^IIoclKMIIz^m) ^ 1, 

by the a priori estimates (|3.9|) and (|3.1ip . Similarly, for the approximately non-atomic term (using 
the L°° estimate of ([311), 

\\wi(s)V • e T/2 v {e T/2 i + z, h e T )\\ LV{jn) = e T \\wi(s)w (e r/2 £ + z i ,e r )\\ LP ( m - ) 

< e T \\w (e T )\\ oo \\w i (s)\\ LP(m) < 1. 

Putting everything together we have 



(r)\\ L 2 {m) < _ + / air- s)- 1 ^ (\\Vw t (s)\\ + K) ds, 

Q-v To) J to 

for some constant K (we have also used that ||iUi(To)||£2( m ) is uniformly bounded). Therefore, for 
some constants Cj, 

a(r - T Q ) l l 2 \\V Wl {T)\\ L , {m) <C l + C 2 f %l T f /2 W2 a(r - ^HV^U)^ 

J T0 a(r- S y/P 
Hence by choosing r G (to, tq + T) for T sufficiently small: 



a(r - r 



1/2 j „ r °( r ) 1/2 J 1 



TQ<T<TQ+T J TO 

and similarly, 



a(r _ s) i/2 a(s _ ro) i/2 Q<T ^ ^ a(r _ s) l/2 a(s) l/2 - 4 > 
sup K j - J ±-L-—ds < 1, 



0<r<f JO a{T-S) 

we have that 1 1 Vi/j^ (''") 1 1 x, 2 (m) 

< 2(Ci + l)a(r-r )~ 1 / 2 for r G (t ,t + T). However, To was arbitrary 
and T was independent of tq so ||Vwi(r)|| i 2( m ) must be uniformly bounded. □ 

By the precompactness just proved, the orbit {wi(r)} T<log , T \ has a non-trivial a-limit set A, 
which we show is invariant under the self-similar PKS (ll.lOp . This is essentially equivalent to 
showing that as r — > — oo the remainder Ri(r) becomes negligible, due to not being localized 
around Z{ . The primary difficulty is the presence of the velocity field coming from the approximately 
non- atomic part wq, which requires some care to properly control (as in |33j). 

Lemma 3. For any i G {1, N}, m > 2 and p G (1, 2), 

lim \\Ri(j)wi(T)\\ Lv{jn) = 0. 
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Proof. Following a similar procedure as [33J, 

||^(r)^(T)(O m || P < UO^Riirn^WMrnmm+i) < IKO^WII^. 

p-2 x ' p-2 

The portion of Ri(r) coming from the other concentrations is relatively easy to handle as we know 
a priori they are localized away from the i-th concentration by the uniform bound on ||z^ 2 (m) • 
Indeed, define q = ^2 an d choose v G (0, 1 — 2/q) (as in [33]). Then by (|2.2p . 

IKfTSfc " (** ~ z*)e- T/2 ^)\\ q < ||<0 _1 « - (** - ^- r/2 >-loolKO^-(OII, 

< e^/ 2 ||^-(-r)|| i2(m) < e^/ 2 . 

Dealing with the approximately non-atomic part is more difficult and is the content of the following 
lemma. 

Lemma 4. For all q G (2, 00], 

hm ||(0- 1 e T/2 Sd(e T ,€e T / 2 +^)k| = Jim t^i ||i)o(t, || L , = 0. 

Proof. We proceed similar to Lemma 4.2 in [33] but with several changes due to the lack of incom- 
pressibility and the unavailability of [19]. Without loss of generality, we can assume Zj = 0. Recall 
that wo(t,x) satisfies the linear advection-diffusion equation (as a mild solution) 

d t w + V • (vw ) = Aw , (3.12) 

where v = B * u, with initial data wo(0) = fiQ. As v satisfies the estimates (|3.9p . the arguments of 
Theorem [T] show also that the linear PDE (I3.12D defines a strongly continuous linear semigroup on 
finite Borel measures which can be expressed in Duhamel integral form 

w (t) = e tA Mo - I >" S ) A V • (v(s)w (s))ds. 
Jo 

The main difficulty posed by the lack of incompressibility at this step is that we no longer have 
the results of Carlen and Loss |19] to provide pointwise estimates on the fundamental solution of 
(13.12p . However, such precise pointwise control is not necessary. 

Since no({0}) = 0, for all 5 > 0, there exists some r > such that ^o(-£>4r) < 5. We may 
define w 1 as the unique solution to (]3. 12|) with initial data /iol_B 4r and w 2 as the solution to A3. 12H 
with initial data ^olR 2 \B 4r an d denote v 1 and v 2 the corresponding velocity fields determined from 
v % = B *w l . By linearity, of course wq = w 1 + w 2 . By properties of A3. 12f) analogous to (i) and (ii) 
of Proposition [6] and Proposition [2] we have for all p G [1, 00] and q G (2, 00], 

1— - 1 

limsupi p\\w {t)\\ p < \\lJ,olB 4r \\M ^ $, 
t\o 

-— - 1 

limsupt 2 i\\v (t)\\ q < \\iM)1b^\\m ^ 
t\o 

Let (j)(x) be a smooth, radially symmetric, non-increasing cut-off function which is one for |x| < 2r 
and zero for \x\ > 4r. Now, further decompose w 2 (t,x) = 4>(x)w 2 (t,x) + (1 — (p(x))w 2 (t,x) := 
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w 3 (t, x) + w 4 (t, x) and v 3 = B *w 3 and v 4 = B* w 4 . One can compute the rate at which mass flows 
into the origin by (using the definition of distribution solution), 

^ w 3 (t, x)dx = — f w 2 (t,x)(fi(x)dx = j w 2 (t, x)A(fr(x) + w 2 (t, x)v(t, x) ■ V<p(x)dx 



dt v ' ' dt 



< u 

< I + IMI. 



Since the RHS is integrable at t = and f w 2 (0, x)4>(x)dx = 0, we have that 

lim / w 3 (t, x)dx = lim / w 2 (t, x)6(x)dx = 0. 
t->o y t->o J 

By interpolation against the a priori L°° bound ||w 2 (t)|| 00 < t , for p £ (l,oo), 

t'-^wm, < t^/^m^iiwsmi < ikwiii -> o. 

Therefore by (|2.ip . for all g E (2, oo], 

limt5 _ 9||u 3 (t)|L = 0. 
It remains to control i> (i,x). For |x| < r we have, 



i i 



*H |^ ( t,x)| < tH / < ^V(*)||r - o. 

7|y|>2r \X-y\ r 

For |a;| > r we have, 

i II 2 2 

|| v 4 (t, a;)i |a!|>r .(i^_)-i || a < Cj)- x £-*\\v\t,x)\\ q -> 0, 

by (|3.9p . Putting the estimates together, we have shown that for all 6 > 0, 

-— - I 2 -) 2 —i 

limsupt 2 9 \\v (t,x)(— — } || < <5, 

which proves the claim by choosing 5 arbitrarily small. □ 

This completes the proof of Lemma [3j □ 

The following is a direct consequence of Lemma [3] and the proof proceeds analogously to Lemma 
6.3 in 1 33 1 so we omit it. 



Lemma 5. A is invariant under the self-similar PKS. 

By standard considerations, A consists only of compact, entire orbits of finite and constant 
self-similar energy Q. However, by Proposition [IJ the G a are the unique functions of this type as 
the self-similar free energy is strictly decreasing on all other sets. Therefore we have proved: 

Lemma 6. A = {G ai (£)} and hence for all m > 2, 

lim ||iOi(r) - G ai \\ L 2f m) = 0. (3.13) 
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We now proceed to step (ii) and construct suitable Wj. As noted above, unlike in [33], f)3.13|) 
does not immediately imply Proposition [3j Indeed, our decomposition is of the form 



N 



N 



i=i 



i=i 



Vt 



Vt 



(3.14) 



but cannot be decoupled into an equation for individual j. Therefore, it is still not obvious how to 
construct the set of Wj for an arbitrary mild solution. Define 



onjWi | log /. 



X Zn 



Vt 



Wi(t,x) - ^G ai 



X Zi 



Vt J' 



(3.15) 



which by (|3. 13|) satisfies linv^.oo ||w)j(r)|| i 2( m ) = in self-similar variables. Note again, that unlike 
the analogous perturbations in [33], Wi ^ Wi. However, it turns out that the perturbations Wi are 
not very far from Wi, so we will produce suitable Wj using a contraction mapping argument around 
Wi in self-similar variables. Precisely, we will construct correctors to Wj, denoted Rj (different from 
above), which will be used to define Wj = Wj + Rj. To do so, we construct solutions {Ri(r, £)}^Li 
to the following system of integral equations 

R i ( T ) = - Wi (r) - 



T ai {r - r')(V • (( Wl + Ri) £(1 - <K£e^2 + Zi _ z . ))v Gj ^ _ {Zj _ Zi)e -r>/l )))dT < 
T ai {r ~ r')(V • C£ —(Wj + RM ~ & ~ z l )e- T '/ 2 me T ' /2 )v G >))dT' 

i+3 

[ T T ai (r - r')(V • (^Gj> G ^ " & " ^ T ' ,2 ))W 

3+i 

J -°° Ul 3+i,3>l 

[ T T ai {r - r')(V • (-G^' 2 ^',^'' 2 + Zi )))dr' 

J-oo a i 

,-z i )e-*/ 2 )))di' 



A' 



T Ql (r - r')(V • ((Ri + Wi) £ <*&(£ - ( 



T at (r - r')(V • ((Ri + W t )e T ' 2 v Q (e T ' 2 + Zi)))dr' , 



(3.16) 



where now Vj = B * (vjj + Rj) and %(t, x) := B * wo. By construction, Wj = Wj + Rj solves (|2,12|) , 
To simplify (|3.16p , first notice that W{ is a mild solution to 

d T Wi + A a Wi = LWi-V- (wm) - V • ((Wi + ^-Gi)(e T / 2 vo(e T ,te T / 2 + Zi ) +% + v 9i fj , 

where v 9i is defined as in (|2,7p and if Vj = B * Wj then we define 

Vi := ^2<XjVj(T,t- (zj - Zi)e~ r/2 ). 

3+i 

Since Wi(t) — > in L 2 (m) as r — >• — oo we can write u), as a solution to the integral equation 



Wi(r) 



tSiUi + (Wi + — Gi) fe r ' /2 z;o(e r ',ee T ' /2 + z t ) + V t + 
on V 



dr'. (3.17) 
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Also write v Ri = B * Ri and 



v r * :=Y,^v Rl (^-(z J -z i )e-^ 2 ) 



Now applying (I3.17P to (I3.16P gives 



Ri(r) = - r T ai {r - r')(V • ((wi + Ri) £(1 - Me" 2 + Zi - zj))v G i (£ - ( Zj - z^' ' 2 )))dr' 
J -°° i+i 

- f T r ai (T-rO(v-(^^(^ + ^o(e-(^-^)e- r ' /2 M^' /2 > G 0)^ 
J -°° i^j i 

- [ T T ai (r - r')(V • {—G i v r *))dT J - f T ai (r - r')(V • (R t (v r > + Vi}))dr' 

J —co J —oo 

T at (r - r')(V • (w iV r >))dT' - f T T ai (r - t')(V ■ (wiv^dr' . 



- f T ai (r - r')(V • (R t e T '/ 2 v (e T ',^ 2 + z t )))dr' . (3.18) 

J —oo 

Lemma 7. For e > sufficiently small and for all T sufficiently small ( as always depending on e ), 
the system of integral equations ()3.18j) has a unique solution {-^(t)}, 1 in C{{— oo, log(T)); L 2 (m)). 
By construction, lim T _ s ._ 00 1 1 -f^i (t) 1 1 z, 2 (m) = 0. 

Proof. The proof is a contraction mapping argument which is essentially the same as the one used 

to prove Propositions [7] and El Using that ||w^(T)||i,2( m ) — >• as r — > — oo most of the arguments 

go forward with little or no changes. The only term which requires special attention is the last 

term in (I3.18D which is being treated here as linear, as vq is being considered as an external field. 

Writing wq as a Duhamel integral involving Swifts) and Wi, it follows from Proposition [6] (ii) and 
i_ i 

()2.ip that limsup t _>o ^ 2 9 ll^o(^)||g ^5 e and therefore should be manageable. Indeed, for p € (1,2) 
using Proposition [5] (iii) 

T ai (r - r')(V • {Riir')^'/ 2 ^',^'/ 2 + z t )))dr'\\ L2(m) < 

e -Hr-r') 



< 



\\R i ( J >)e'/ 2 v (e r ',Ze'' 2 + Zi)\\ L P{m) dr' 

j a(r — t') l /p 

e -v{r-r') 

a ( T _ T ,y/p \\ eT ' /2 Me T \& T ' /2 + Zi)\\2 p /( 2 -p)\\M T ')\\LHrn)dT. 



However, (writing t = e T ), 
By Proposition [U 



e T ' /2 v (e T ' ,te T ' /2 + z i )\\ 2p/(2 _ p) =t^ V \\v (t,x)\\ 2p/{2 _ p) . 



1 2-p 

limsupt 2 2 p \\vo(t,x)\\ 2p /(2-p) i$ WvoWpp < e. 
Hence, for r < log(i) for t sufficiently small, we have, 

II f r Qi (r-r')(V-(i2 i (r') e -'/2 U0 (^72 + 2; . )e r' )))dr /|| i . 2(m) < eM[i2i](e -). 

J — oo 

The rest of the terms can be handled using the same techniques as the proofs of Propositions [7] and 
Accordingly, the lemma follows from the contraction mapping theorem. The decay as r — > — oo 
follows directly from (|3.18p . □ 
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There is one last remaining detail: due to the nonlinear terms in (|3.18p it is not immediately a 
priori obvious that one re-constitutes the original ^2iWi(t,x) by summing the Wj. From 

El ( X — Zi\ 1 / X — Zi\ 1 / X — Zi \ , . 

:1 i Ga > {~7r) + ~t Wi l logt ' ~^r) = h 1 v ~7r ) + wdt ' x) ' 

it suffices to prove the following lemma. 
Lemma 8. For all t € (0,1% 

N 



i=l 

Proof. By definition, 

N 1 fx-zA N 1 / x-Zi\ N 

^2wi(t,x) + ~G a ( — -j^- J = ^2 J R i ( lo g^ — ~7t~ J +Wi(t,x) := F(t,x) + 

i=l i=l i=l 

where F(t,x) := Y^i=i ( x ~ Zi)^ 1 ! 2 ). By construction, f(t, x) = F(t, x) + X)i=i w i(t> x ) 

is a mild solution to 

d t f + V-(f(v + v F )) = Af, 

with initial data /(0) = Y2i=i where v(t, x) = B*u, the original velocity field, and v F = B*F, 
the potential defect introduced by F. By definition, w i solves the advection-diffusion equation 
with the same initial data and velocity field v(t, x) and therefore F is a mild solution of 

N 

d t F + V • {Fv F ) = AF - V • (Fu) - ^ V • 

i=i 



with initial data -F(O) = 0. This is equivalent to 

N 

V • (F( S )v F (s)) - V • (F(s)v( S )) Ms)v F (s)) 



F(t) 



* (ts)A 



t=l 



However we know a priori that lim^o * 1//4 ||^(*)ll4/3 = (since linv-^-oo [ | [ I i, 2 (-m,) = 0)- Using 
the above integral equation, (|1.6p and the a priori estimates on v and W{, it follows that for t 
sufficiently small 

t 1/4 ||F(t)|| 4/3 < \ sup S 1 /4|| F(s) || 4/3) 

which implies that F(t) = on some small time interval. Iterating implies that this holds for all 
t > such that the a priori estimates on v and u>j hold. □ 

Since Y^iLl r i( t ^ x ) = °> 

JV N 



X> = E (^) + ( lo ^ x -jf) 



i=l i=l 

AT 



1 /"» f x ~ z i\ . a i ~ (i X- Zi 

Therefore, the above construction yields a proper decomposition of the original mild solution u(i, x) 
Moreover, Wj satisfy (|2. 12|) . This completes the proof of Proposition [3l 
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3.4 Final Step of Theorem [2] 

Proof. (Theorem [2]) By Propositions [7] and [U the contraction mapping theorem implies that there 
exists a unique solution to (|2.10p and (|2.12p which satisfies M[w](t) < e for t small. By construction, 
we may assemble a mild solution u(t,x) to PKS via (|2.3|) . By Theorem [H this mild solution also 
satisfies the L°° estimate \\u(t, x)||oo < i -1 for short time and due to the regularization effects of 
the linear diffusion, after t > the solution is classical for as long as it exists. Proposition [3] implies 
that any other mild solution can also be written as a solution to (|2.10p and (]2.12p and, by the 
uniqueness implied by the contraction mapping theorem, must agree with the constructed solution. 
This completes the proof of Theorem [2J □ 

4 Lipschitz Dependence on Initial Data 

Although the contraction mapping theorem generally provides Lipschitz dependence on initial data 
for free, the above argument does not. This is because all of the linear propagators used in the 
proof of Theorem [2] depend on the initial data itself and the decomposition used to construct the 
solution. However, combining our decomposition (12. 3p and the contraction mapping arguments of 
Section ^3.1l with the main ideas of Section 5.3 of (33] , we are able to refine the results of Gallagher 
and Gallay to obtain Lipschitz dependence on initial data. As the method applies equally well to 
the 2D Navier-Stokes equations in vorticity form (jl.2p . we also prove that the solution map of the 
Navier-Stokes equations is locally Lipschitz. 

4.1 Proof of Theorem [3] 

Proof. Let e > be a small, fixed constant independent of 5 to be chosen later. In particular, we 
can require 6 to be small with respect to e. If 5 < e then all the atoms of mass bigger than e in /i 1 
and ii 2 must be in the same location (this is the utility of the total variation norm as opposed to a 
weaker norm). Therefore, as in |33| . we may decompose the initial measures as 

N 

such that ll/^ollpp < e and 

N 

ii^ 1 - v 2 \\m = H/4 - i4\\m + ^ \ a } - a )\ 

3=1 

Let u , u 2 be the unique mild solutions of (jl.ip in \T associated to each of the respective initial 
measures constructed in Theorem [2j By Proposition [3l we may decompose u l each into respective 
large atomic pieces and smaller perturbations as in (|2.3[) : 

u l (t, x) = w l (t, x ) + J2 ] G 4 (^tt) + ^ ( log *' ^7r) ' (4,1) 

where the perturbations w\ satisfy the corresponding integral equations (|2.10p and (|2.12p . Analo- 
gously to [33j we use the following quantity as the norm to measure the difference between the two 
solutions, 

A(t) := sup \\wq(s) -«>jj(s)||i +M[w l -w 2 ](t), 
se{o,t) 
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where for future convenience we define 



A (t):= sup ||^(s)-^(s)|| 1 + K MoK-^](t). 

se(o,t) 

By Proposition [1] and the decomposition (|4.ip , Theorem [3] is equivalent to: there exists a T > and 
Cl < oo (independent of 5 and // 2 ), such that A(T) < Cx<5. 

The decompositions of define different linear propagators which are centered around the 
same points but have different masses in the concentrations, denoted S l N (t,s) and 7^i(r). The 
difference between these linear propagators must be controlled, hence in order to continue we need 
the equivalents of Proposition 5.5 and 5.6 of [33]. The proofs are a straightforward contraction 
mapping argument which follows by writing one linear propagator as a Duhamel integral equation 
involving the other and using (iv) and (v) of Proposition [1] to estimate the terms in the integral. We 
omit the details of the contraction mapping, but the proofs of Proposition [T] (iv) and (v), which are 
trivial for NSE but not for PKS, can be found in Appendix fjBl For Sw(t, s) we need the following. 

Proposition 9. There exists some to sufficiently small (depending on d and ||/i ||ppj such that the 
following holds (with implicit constants independent of 5 and e): 



(i) For all p G [l,oo] and v E Ai(M?) we have 
|| {S 1 N (t,s)-S 2 N (t,s))u 



I LP < 



(t - s y- i /p 



\v\\m, 0<S<t<S + t(). 



(4.2) 



(ii) There exists some Ao £ (0, 1/2) independent of e such that the following holds: for all 7 > 
sufficiently small and for all f E L 1 (IR 2 ) ; 



(Sh(t, s) - S%(t, s)) V/|| p < {f _ s) 6 3/2 _ 1/p ' 1 



7+1/2-A0 



1, 0<s<t<s + t . (4.3) 



Similarly we need the analogous estimate for T a . 

Proposition 10. Fix a E (0,87r), m > 2. Then for some v G (0,1/2) (which depends on a), all 
5 6 [1,2], all (3 sufficiently small (in absolute value), and all f G L q (m), 



(T a+ p(r)-r a (r))Vf\\ L2{m) < a , q |0| 



where the implicit constant is independent of f3. 



a(r)V« 



L«(m)> 



r > 0, 



(4.4) 



Armed with these estimates we may now use arguments similar to those in Section ^3.11 to 
estimate the norm A(t) using the integral equations satisfied by w\. Consider first, 

wl(t) - w 2 (t) = (S* (t,0) - S 2 N (t,0)) ill + S 2 N (t,0)(»o ~ Mo) 



(S 1 N (t,s)-S 2 N (t,s))V 



N 



ds 



S l N (t,s)V 



N 



N 



3=1 



3=1 



ds. 
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An argument using the known a priori estimates on u l proves that (for to sufficiently small so that 
Propositions [6] and [9] hold) 



A (t) < SK^l + M[w 1 ](t) 2 ) + tf 2 (M[w 1 ](t) + M[* 2 ](t)) A(t), < t < t 0j 



(4.5) 



for constants K\, K2 independent of e. We now turn to the estimates for w^, i 6 {1, iV}. In what 

k an< ^ 



follows write Gj^ = — and 



4 



(r,0 :=^a^.(r,e-(^-^)e- r/2 ), 



-t/21 



(4.6) 



Consider the first term, 



^ T a i(r - r')V • (^M) " T a 2(r - r')V • (^M) 



r a i-r a? )(r-T0v 



5^ I dr' 

tts 



TAt-t') V- 



07 



G, 



r Q i-r a? )(T-T0v- 

G„i 



^L^ 1 ) r/r' 



r T ai (r-r')V 
J — 00 



07 



tP i — v 



,91 



dr' 



G n 2 G„i 



07 



O ; 



dr'. 



By (1X211 and gaj), 



e -^(r-r') G a i , 
\Gi,i(r)\\ L 2 {m) < I <5 a(r _ r / ) i/ 2 W^t v9i \\l?{m)dr' 



at 



+ 



+ 



(r - r') 1 / 2 " a, 1 



ySi _ v 9i) || L2(m) dr / 



-v{t-t') 



a(r - r') 1 / 2 



G. 



G., 



07 



- I ^ llL2( m ) 



dr'. 



Notice that while the third term is always zero for the Navier-Stokes equations, it is non-zero for 
PKS due to the nonlinear nature of the self-similar solutions. The first term can be estimated as 
we estimated Fn in Lemma [TJ which gives the following with an implicit constant independent of 
e and Oj, 

e-Kr-TO G a i 1 , .„ 
U ^rV^\\ L2{m) dr' < bfl\ 



a(r — r') 1 / 2 a 

The second term can be estimated in a similar fashion but now using Proposition Q] and Proposition 
[2] to deduce the following with an implicit constant independent of e and aj, 

e ~ U{T ~ T ' ] " ' ^ ^ '\ LHm) dr' < W 2 \a) - a]\ < e^8. 



i(r-r') 1 / 2 " aj 
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Finally the third term can be estimated similarly with an implicit constant independent of e and 
Proposition [Q 



— V y < || L 2( m) 



a(r — r') 1 / 2 \ a 



at 



~\\L 2 (m+l) 



p t/2 „r/2 

< e — | a ? _ Q i| < f_,5 



Hence putting the three estimates together we have, 

l|Gj,i||L2( m ) < Se r/2 ^1 + . 

Dealing with the other terms does not really present any new real challenges, so we will include less 
details. The term Gj 2 is done in a similar way: first re-write as 



G , 



Gi,2(j) =J2 I T a] {r - r')V • dr' + £ I T a ,{r - r')V • I ^)v)\ dr' 



G a 2 



+ E f T^r - r')V • {% [a)v) - a 2 v 2 )) dr> 



dr'. 



By estimates analogous to those used for F^i in Proposition [7] combined with Propositions Q] and 

EDI 

\\Gi,2{r)\\ L 2 {m) < 6e / / 2 ~ 1 / q M[w 1 ](t) + ^e^-^M^t) + A(t)e^ 2 -^, 

at- 
1 

where 7 and q are as in the proof of Lemma [1] and the implicit constant is independent of e and 5. 
Similarly, we re-write the third term as 

Gifl = (T a} - T al ) (r - r')V • (^7%V, £ e *7* + ^ ^ 



+ / r a? (r-r')V 

j —00 

+ f T 2 (t-t')V 



1 2 



^^' /2 ((^--o 2 )(e^,^'/ 2 + ^)) 



dr' 
dr'. 



Using ideas from Lemma Q] we deduce that there exists some K3 independent of e and i such that, 
\\G^(r)\\ L 2 {m) < K z 8M Q [wl\{t) + K 3 ^M [wl](t) + K 3 M [v% - w 2 ](t). 
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We may continue in a similar manner to handle the remaining terms and eventually prove an 
inequality of the form for some r\ > using also that M[w l ](t) is uniformly bounded, 

IK- - WihHm) <K A \1 + + -2 J 5 + K 3 M [wl - w 2 ](t) + K 5 (6 + e" T + M[u) ](i) + M[* 2 ](i)) A(t) 

< iT 4 f 1 + e" T + t4^2 1 ^ + ^A(*) + ^5 + + M[w l ](t) + M[w\t)) A(t), 
1 mi a - I Kq 



3 



(4.7) 



where K4, K3 and K§ are all independent of e and <5. Choosing e > small depending only on 
absolute constants, Kq > max(l,4K3) and t < to for to sufficiently small depending on e. Note, 
none of these choices depend on S as long as 5 < e. Moreover, we may T small such that M[tu'] < e 
by the work of Section 52 Taking the supremum over i, from (|4.5p and (|4.7p we deduce that 

A(t)<5, < i < i , 

which as remarked above, implies Theorem [3j □ 
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A Appendix: Linear Estimates 

A.l Proof of Proposition [5] 
A. 1.1 Proof of (i) 

Fix a £ (0,87r). Consider / which solves d T f 
Duhamel's formula this is 



= Lf — A a f with initial data /q. Written with 



f{r) = 5(r)/ - / S(t- s)A a f(s)ds. 
Jo 

A straightforward contraction mapping argument similar to the others employed in this work implies 
that 7^(t) = e T ( L --^a) d e fi ne s a strongly continuous semi-group on L 2 (m). However, an additional 
analysis must be done to ensure (12.180 holds independently of r. By linearity, it suffices to consider 
nonnegative /o with ||/o||L 2 (m) = 1- As above, write f(r) = T a (r)fo- As L 2 (m) <— > L , /(r) € L . 
Furthermore, 7^(r) preserves non-negativity and is in divergence form, hence ||/(t)||i = ||/o||i f° r 
all r > 0. 
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< 
< 



We first show that ||/(t)||2 < II/0II2 independent of r. Indeed, 

= - j Wf?dt + \ j \f\ 2 di + ^ j G a \f\ 2 <%- j fVG a -Vcd£ 

Nf\ 2 dt + c\\f\\ 2 2 + c\\f\\ 4/3 \\vc\u 
\vf\ 2 dt + c\\f\\ 2 2 + c\\f\\l /3 
<-J iv/i^ + qi/nl + qi/iixii/ib 

<-J Nf\ 2 dt + C\\f\\ 2 2 + C\\fo\\l 

for some constants C which depend on a but whose precise values are not very relevant. Similar 
to the proof of Theorem [1] (i) , we apply the Gagliardo-Nirenberg inequality \\f\\^ < IIV/H2II/H1 and 
get 



~ j \f\ 2 dt<-c^ + c\\f\\ 2 + c\\f \\l 



Using that for all K > 0, WfWl < jj^ll/lli + C(-fQII/l|i the above differential inequality implies that 
||/(t)||2 is uniformly bounded by a constant. 

Now we may use the uniform bound on ||/(t)||2 to control || |£| m /H2, similar to what is done in 
Theorem 3.1 in [35]. Computing as there, 

lit S l^l 2m ^ = / \i\ 2m f{Lf-Kf)dt 

= - I l?| 2m |V/| 2 ^-m J \i\ 2m f 2 di + 2m 2 j \tf m -*fdi-j |e| 2m /A Q M. 
The latter term expands to the following after a short computation (using that —Ac = /) 

-j \i\ 2m Kfdt = \J \i\ 2m G a f 2 di- f |e| 2m /(VG a ■ Vc)de + m j \t\ 2m - 2 i-vc a f 2 dt. 
By the rapid decay of G a the first term is uniformly bounded via 



. \zrG a m<\\f\\i 2 . 

Using Holder's inequality and the L 4 estimate (|2.1j) we control the second term 
- J |e| 2m /VG Q -Vc^< |||e| 2m VG Q || 4 ||/|| 2 ||V C ||4 

<l|/l|2||/||4/3<ll/||2||/|k 2 (m)- 

For all 8 > there exists a constant C$ > such that the following three inequalities all hold (using 
the spatial decay of Vc a ), 

\\fh\\fh*(m) <*I|/Il| 2(m )+Q||/Ill, 
2m 2 J \i\ 2m ~ 2 f 2 di<5j \i\ 2m fd£ + C s J f<%, 

m j \t\ 2m - 2 i-Vc a f 2 dt<8 j \H\ 2m fdZ + C s J f 2 dt 
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Putting all of the estimates together we have, 

\tt J \tf m f 2dT ^-J ie| 2m |V/| 2 ^-m J \H\ 2m fdH + 35 J \tf m fdH + C J M. 

Therefore, since ||/||2 is uniformly bounded, for 5 chosen sufficiently small this inequality implies 
that || |£| 2m /H2 is also uniformly bounded. Since 7^(t) is a linear operator the bound (|2.18p must 
hold. 

A contraction mapping argument similar to what is done in Lemma 2.1 in |36j and Lemma [2] 
shows that the uniform bound (|2.18j) implies the regularization estimate (|2.19p . Moreover, in both 
cases, one can verify that Proposition [T] (iv) implies the implicit constants in (I2.18D and (|2.19p only 
depend on K < 8tt for a < K. Note that T a can be treated as a small perturbation of S(t) for a 
small. 



A.1.2 Proof of (ii) 

J. Campos and J. Dolbeault have recently proved the following spectral gap. See also Appendix O 
for an independent proof. 

Theorem 4 (Spectral Gap (J. Campos and J. Dolbeault p2])). F° r a ^ a G (0, 871"), there exists 
some K a S (0, 1/2] such that any eigenvalue A of L — A a in L^G" 1 ^) satisfies 

Re(A) < -K a , 

for some K a > which is uniformly bounded for a small. 

We extend this result to the polynomial- weighted spaces in a way analogous to [36] : 
Proposition 11. Fix m > 1. Then any eigenvalue A of L — A a in L^m) satisfies 

(2 — ffi 
—K a , — — 

Proof. As in section 4 of |36| . we use an ODE argument. The idea is to show that any eigenfunction 
of L — A a in Lq(7tt.) actually lies in Lq(G~ d£) and hence we may apply TheoremHl The approach of 
Gallay and Wayne [36] is to re-write the linear operator in radial variables and reduce the question 
to a statement about the asymptotic behavior of an ODE, for which classical results of Coddington 
and Levinson may be applied [25]. Hence the core of the argument is the following lemma. 

Lemma 9. Let f € L^{m) satisfy Lf — A a f = [if for Re(/it) > Then there exists 7 > such 

that 

1/(01 <(i + iei 2 ) 7 e- |f|2/4 , erf. 



Proof. As in [36J, we decompose the eigenfunction with a Fourier transform in the angular variables. 
Each mode decouples due to the radial symmetry of the coefficients. Define —Ac = / and write 
/ = Yl^=-oo fn{r)e tne and c = ^^-00 c n (r)e md . Written like this we have 

00 r 1 / 2 \ 
Lf-A a f= £ e me -(r/;)'+(^-4(r))/;+M+2G a (r)-^J/ n -GUrK(r) 

n=— 00 L \ / . 

1 , , n 2 

fn = i rc n) H 2 Cn ' 



38 



Due to this decoupling and linearity, we may assume without loss of generality that the eigenfunction 
is supported in only one mode, n, and writing f(r) := f n {f) and c(r) := c n {r) we have 

- r {rf )' + g - c») f+(l-n + 2G a {r) - ^ / - G' a (r)c'(r) = 

1 n 2 
— (rc')' + -5-c = /. 

As in [36] we will re-write this ODE in a form amenable to a classical result on ODEs (Theorem 
8.1, pg 92 [25], and following Gallay and Wayne, we introduce a new set of variables 



~.2 



t = j, f(r) = g(t). 



In these new variables we have 



/ 1 d a {2y/t) \ /1-A* m G" a (2v^)c'(2y^) 
5 (*) + I 1 + ^ J 9(t)+ I — a(t) I = , (A.l) 

where 

n 2 2G a {2yft) 

a{t) = tf 1 • 

We first analyze the linearly independent solutions of the homogeneous equation, 

g"(t) + ~9'it) + (±=± - o(t)) git) = 0. (A.2) 

The primary way this ODE differs from the corresponding one for the Navier-Stokes equation is the 
term involving c' a , which is only present for PKS. Note that 

lim 2Vtc' a (2Vt) =-■£-. 

Now, define x(t) = (g(t),g'(t)) and re-write (|A.2|) as the system 

x'(t) = (A + V(t)+R(t))x(t), 

with 

A =(S -i)< ^)=(_L fl(4) =U 2)- 

Hence by Theorem 8.1, pg 92 of [25] we can get information about the decay of solutions by analyzing 
the eigenvalues of A + V(t), given by 



1 1 , c' a (2Vt) 1 // 1 c^/i)\ 4-4 M 

Hence lim^oo A + (t) = and lim^oo A_(i) = —1, the eigenvalues of A In order to use the result 
of Coddington and Levinson, we only need to compute the terms which are unbounded in T in the 
time- integrals 

T 

\±(t)dt, 
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as the integrable terms can be absorbed into the constants. Using a Taylor expansion of the square 
root, dropping the terms which are integrable as t — > oo we get 



Ah 



I, M i 

2t 2^1 4 
1 c> a (2Vi) 1 
2* 2^/t 4 



2 4^ 2<4(2y/f) 
2 4 M 2c4(2Vt) 



Simplifying, 



A_ 



1 _/x + c' Q (2^) 



■1 



/i a 
7 ~ Zvr? 



The theorem of Coddington and Levinson then implies that there exists two linearly independent 
solutions to (jA.2|h 4>i(t), 4>2(t), such that 



lim t 1 -^ 



t— >oo 



/ ' t—>oo 



lim t^ + ^e 



lit) 



1 



Returning now to inhomogeneous ODE (jA.ip we use the variation of constants formula, as in |36j . 
Note that the inhomogeneity here is 



b(t) 



G' a (2Vt)c'(2Vt) 



t 



Writing the solution g(t) with the variation of constants formula gives 

g(t) = A(t)Mt) + B(t)Mt), 

with 

rt 

\-l 



A(t) = A x - j (W(s))- L b(s)Ms)ds, B(t) = B x + J (W( S )y L b(s)Ms)ds, 

where the Wronskian W(t) = 0i(t)<^(i) - <fo (i) « -i -1- ^ - * as t — )• oo. Since / has average 
zero, 

| c '( r )| ^ — > r ^ oo. 



Hence by Proposition [U 
which implies that 
As t — >■ oo, 



l^)l<*" 2 e"*, i^oo, 
WC*)- 1 |6(t)| < t^- 1 , t->oo. 



w{tr i \b{t)\\ut)\<t~ l -^~\ 

which is integrable, so A(t) — > A^ for some constant. On the other hand as t — > oo, 

W{t)- 1 \b{t)\\4> l {t)\<t^+^ 



which implies 



m)\ <(!+*)* 
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By the asymptotic behavior of <j>\, for large t we have 

Returning to the original variables shows that ^ would violate f £ L 2 (to) and hence must be 
zero. By the asymptotic decay of fait) together with the polynomial bound on B(t) we have the 
claimed exponential localization. □ 

The lemma shows that any eigenfunction / £ L^{m) with eigenvalue Re// > (1 — m)/2 must in 
fact be in L 2 {G~ d£), and the result follows by Theorem HI □ 

The next step is to prove the decay estimate for the time-dependent linear evolution equation, 
for which we follow the arguments of section 4.2 in [36] . In what follows, for any linear operator 
A, denote a ess (A) the essential spectrum, defined as the set of A 6 cr(A) such that XI — A is not 
Fredholm (in some textbooks this also include A such that XI — A is Fredhom with i(XI — A) ^ 0). 
Recall that compact perturbations do not change the essential spectrum [32]. The essential spectrum 
of 5(r) in Lg(m) is known to be {A G C : | A| < e^ 1 "" 1 )/ 2 } (L emma 4.11 [36]). The following lemma 
shows that T a {j) is a compact perturbation of S(t). This ensures that the essential spectrum of 
T a {r) is the same as S(t), and so for m chosen large, only the spectral gap given in Proposition [TT1 
is relevant for estimating the time decay. 

Lemma 10. Let m > 1. The linear operator K{t) = T a (r) — S(t) is compact in L 2 {m) for all 
t > 0. 

Proof. Let /o £ L 2 (m) and write 

K(r)f Q = - f T S(r-s)A a f(s)ds, 
Jo 

where /(r) := T a (r)f . Then by $ZM), 



\\ K (i-)fa\\L2( m +i) < / \\S(t - s)A a f(s)\\ L 2 {m+1) ds 



e -§<V-s) , f \ 

1 1 a f T _ s )i/2 [)\ v /( s )iU 2 (m+i) + l|ga^(g)[|^(m+i)j da, 



where v* = B * f. To control the first term we use the spatial decay of v Ga , 

ll^/llx^+l) < ll« Ga (0<Olk«ll/(s)llL»H < ll/OOIIx 2 ^) 

The second term we use the L estimate (|2. 1 1) . 

Ik G'a 11x2(771+1) < Ik llx^l^allx^m+l) ~ ll/llx 4 /3 ^ 1 1 / 1 1 X 2 (na) • 

Hence, 

/•t e -i(r-s) 

ll^( r )/o||x2(m+i) ^ / -; rT^||/(*)||ia( TO )da. 

However, since 7^(t) is bounded on L 2 (m) we have 

l|£a( r )/o||x 2 (m+l) < r ||/o||L2( m ). 

The estimate (|2.19p similarly implies 

||^a(7")/o||frl(m) ~ ll/o||x 2 (-m)- 

Compactness follows from the Rellich-Khondrashov embedding theorem. □ 
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Following exactly the argument of Lemma 4.11 in [36] we have that the essential spectrum in 
Lg( m ) of 7~(t) is given by 

^(r Q (r)) = {AGC:|A|<e^- m )/ 2 }. 

Combined with the spectral gap estimate in Propositioning we obtain the equivalent of Proposition 
4.12 in |36| with the same proof: 

Proposition 12. Let v € (0, K a ) and /o £ L/^(rn) for some m > 1 + 2v. Then, 

||r a (r)/ || m < v e~ UT \\fo\\ m . 

A. 1.3 Proof of (Hi) 

The proof of (hi) is analogous to Proposition 4.6(iii) in [33], and one can easily check that the proof 
contained therein can be carried out in our case without any significant changes. However, notice 
that the spectral gap estimate deduced in (ii) is vitally important. 

A. 2 Proof of Proposition [6] 

Parts (i) and (ii) can be proved with methods similar to those used in Theorem [TJ However, by 
linearity there is no size restriction on the initial data and the dependence on the norm of the initial 
data is of course linear. More precisely, we construct the solution u(t) := Sw(t,s)v for short time 
using the contraction mapping theorem and the Duhamel formula 

S N (t, s)u = u(t) = e^ A u - J* e (*-*') A V • {v{t')u(t'))dt' , (A.3) 
where we have denoted 

Indeed, using that \\v(t)\\ 4 < t~ 1/4 and (fL6|) . 

(t - s) 1/4 ||u(i)|| 4/3 < Mm + (t- s) 1 / 4 J\t - tT^Ht'Ht^hdt' 

< \\u\\ M + (« - ( sup (f - s)^\\u(t>)\u /3 ) [\t - o- 3/4 (*T 1/4 (i' - sY^dt'. 

\t'e(s,t) J Js 

The latter term is linear however, t can be chosen sufficiently small to make the coefficient small, 
and a contraction mapping proves that one can construct a solution to Sjsf(t,s)u for s > which 
satisfies (I2.22D for p = 4/3. A bootstrap similar to Theorem [JJ (i) shows that this implies the other 
estimates contained in (|2.22[) . Classical theory can then be used for t > s to extend the solution for 
longer times. Inequality (12.23H follows from the same argument and the corresponding property of 
the heat kernel f)1.7|) . 

The remainder of the section will be devoted to the proof of (hi) which is more difficult. We 
proceed similar to Appendix 6.2 of [33], however, the core of our argument needs to use some theory 
on the spectral gap of linear Fokker-Planck operators, due to the gradient nature of the nonlocal 
velocity law. As in Appendix 6.2 of [33], it is easily seen by applying (|2.22p that it suffices to prove 
(|2.24p for p = 1, so we concentrate on this. 
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A. 2.1 Part (Hi) with N = 1: One Concentration 

As in |33j . we proceed first with the case of a single concentration, which requires some known 
spectral gap properties of Fokker-Planck operators. We study, for fixed a G (0, 8tt), the PDE 

fll „ + V.(-Lv C<> (^)»)=Aw. 

In self-similar coordinates, £ = xt -1//2 , r = logi, tw(t,x) = /(r, £), we may re-write this as 

r / + V • (fVc a ) = Lf. (A.4) 

Denote by S\(t) the linear propagator to (|A.4p . The PDE ()A.4j) . for r > 0, is a linear Fokker-Planck 
operator with confining potential 

A(0 :=i|e| 2 -c a (0, 

and it is known that operators of this form admit a spectral gap in the appropriate weighted space, 
in this case L 2 {G~ l d£). The easiest way to see this is to re-write (IA.4D as a self-adjoint Schrodinger 

— 1/2 

operator for the variable z = fG a and apply a classical result that such operators has a pure 
discrete spectrum under certain conditions which are satisfied here (Theorem XIII. 67 |56j). See e.g. 
m for a more detailed explanation. The existence of a spectral gap implies that there exists some 
X a > such that for all f G L 2 (G~ l dCj with / f Q d£ = 0, 

\\m\\ LHG -i dx) < e- AQT ||/o|| L2(G - ld0 , for > 0. (A.5) 

We will assume that X a G (0, 1/2). We now show that (jA.5j) ultimately implies 

l|gi(r)V g /o[|i< 6 ^/ ll/olli, (A.6) 
for any 7 > sufficiently small. When we transform back into physical coordinates, (|A.6|) becomes 



\Si(t,s)V x w\\i < 



^-l 7+1/2- A a 

Iklll) 



1 

noting carefully that we are taking r = logt — logs. To prove (|A.6j) we proceed analogously to [33] , 
We define the Banach space X, 



X = {fe L^M 2 ) : / = d xx g x + 8 X2 g 2 , 9l ,g 2 e L\R 2 )} , 
equipped with the norm 

\\f\\x = ll/lli + inf {|| 5 i||i + \\g2h : / = d Xl9l + d X2 g 2 } . 
As in [33], consider the auxiliary equation for g = (gi, g 2 ) 

d T g + Vc Q V ■g=\L-^\g, (A.7) 

and denote the associated linear propagator by T\(t), which is related to Si via 

V-(T 1 (r)g) = S 1 (r)V-g, (A.8) 

for g £ (if 1 (M 2 )) 2 . A contraction mapping argument similar to Lemma 6.4 in [33] shows the 
following. 
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Lemma 11. T\(t) defines a strongly continuous semigroup on L 1 (M 2 ) 2 and there exists some tq > 
such that for all g G L^K 2 ) 2 , 

WTxi^gh < \\g\U, HVTi^lli < — ^ Nli, r€ (0,7ft]. (A.9) 

By (|A.8p . Lemma [TT1 implies that for r G (0, To), 

l|Si(r)V/||x<— ^ll/lk. (A.10) 

To prove ()A.6|) . it suffices to verify the spectral gap- type estimate 

\\Si(r)f\\x<e^ + ^\\f\\x, (A.fl) 
since (jA.lf | and ()A,10p together imply for r > To, 

ll^iWV/IU < e^+^-^WS^fWx 

e (-A a +7)r 

~ / vT/9 ll/lll' 
a(T) L ' z 

To prove (jA.lip we follow a procedure similar to what is carried out above to prove (|2.20p . By 
writing the solution /(t) = 5i(r)/o as the integral equation 



f(r) = 5(r)/ - f 5(r - r')V • (/(rOVc^dr', 
J o 



similar to Lemma [TUl 5*i(r) can be shown to be a compact perturbation of S(r), which has spectral 
radius e _T / 2 in X [33]- It then suffices to show that all the eigenvalues of the Fokker-Planck operator 
Lf — V • (/Vc a ) have real part less than — X a (recall we are assuming without loss of generality 
that X a G (0, 1/2)). Suppose that there exists some w G X with 

Lw - V • (wVc a ) = fiw (A.12) 

for some ji G C with Re/x > — A Q . As in the proof of (]2.20p . by the radial symmetry we may can 
assume w(r cos 9, r sin#) = f(r)e ind for some n£Z and re-write (|A.12p as an ODE for f(r). Similar 
to the argument in Lemma [9] (easier as the ODE can be treated as homogeneous) , w G X satisfying 
(lA~T2l> implies 

f(r) ps Kr 2 ^- 2 + 0(r p e~ r2/4 ), r -> oo, 

for some p > and if G C However since w £ X and Re/u > — X a > —1/2 it follows that K = 0. 
Therefore, u; G Lq(G~ «£), but this would contradict the spectral gap (IA.5P for the Fokker-Planck 
operator. Hence, necessarily the spectrum of Lf — V • (/Vc a ) in X must be contained in the set 
{A G C : Re A < — X a }. This in turn implies (jA.lip which completes the proof of (IA.6P and hence 
(|2.24p in the special case N = 1. 

For the next section we also need the following perturbation lemma which shows that X a ~ 1/2 
for a small. This is important to ensure that Ao in ()2.24p can be taken independent of e. 

Lemma 12. For all 5 > sufficiently small, there exists as such that if a < as then S\(t) satisfies 

\\Si(r)f \\x <s e ( - 1/2+<5)T ||/o|U, 
where the implicit constant is independent of a. 
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Proof. Let fo £ X and write /(r) = <Si(r)/o in the Duhamel integral form 

fir) = 5(r)/ - T 5(r - r')V • (/(rOVc^dr'. 

JO 

Using the known spectral gap of S{t) in X, and Proposition [1] (iv), 

e T/2 - 5T \\f(r)\\x < e T / 2 - 5T \\S(r)f \\ x + e T / 2 " fr f \\S(r - r')V • (f(r')Vc a )\\ x dr' 

Jo 

<e- 5r \\M\x + e^ 2 - 5T [ T e-^Wfir'^c^hdr' 

Jo 

<e- 5 / 2 ||/ ||x+af snp e^ 2 ^!/ (r')|| x ) ~ (e » - l) . 

\r'e(0,r) J \ J 

Taking the supremum in r of both sides the lemma follows by choosing a < 5. □ 
A. 2. 2 Part (Hi) with N > 1: Multiple Concentrations 

As in [33], to extend to multiple concentrations, we use the intuition that if t/d 2 is small, then 
separated concentrations should basically decouple. Introduce a nonnegative, smooth cutoff x( x ) 
which is one for \x\ < 1/2 and zero for |x| > 3/4. The localizations around the corresponding 
concentrations are Xi( x ) = X (^ir 1 )- Define the opposite localization to be Xo( x ) = 1 — YliLi Xi( x )- 
Proceeding as in [33J, if f(t,x) is a solution of (|2.9p . define fi(t,x) = Xi(x)f(t,x) and note that for 
ie{l,..,N}, 

dfi „ / 1 r ( x — Zi 

with 



x) = ±=vf (^-jf\ Xi( x ) + 2V Xi (x), 



Qi(t, X ) = J2 ^r t v f (^Tf 1 ) ' Vxi ^ + A ^ x )- 

For i = we have (with the remainders Qq, Ro defined similarly), 

^ = A/ + Qo/-V-(i?o/). 
Note that the a priori estimates on /, Proposition Q] and the definition of the cutoffs imply 

N N 
i=0 i=0 

for all t > 0. Denoting S l (t,s) the linear propagator associated with the concentration centered at 
£j we have, for t > s > 0, 

/<(t) = S\t, s)fi(a) + / S^t, t') (Qi(t')f(t') - V • (Ri(t')f(t'))) dt'. (A.13) 
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By the previous section, we know that each Si(t,s) satisfies 

7 +l/2-A aj 

\\ w \\i- 

^From here we may proceed as in Proposition 4.3 of [33J to finish the proof of the Proposition [6] 
with Ao := mini<j<jv(A Q , i ). Lemma [T2l shows that Ao E (0,1/2) uniformly in N, since only large 
values of a can have a relevant effect on A a . ,/From the proof we see that we need to choose to in 
Proposition [6] such that 

< mm(l,K) , 

where K is some constant which is independent of e, ./V and d. 



\Si(t, s)Vw\ 



< 



(t 



a/2 



B Appendix: Properties of Self-Similar Solutions 

B.l Sketch of Proposition [I] 
B.l.l Part (i) 

Biler et. al. proved the uniqueness of G a in radial variables in [7J by an ODE argument. The 
existence of G a had been studied previously in, for example, [52] . ^From the definition of G a , 
for a given mass a € (0, 8ir), G a satisfies the equation 

e c«(0-|£| 2 /4 
G n (£) = a 5 = — A Cr.. 

^From this definition, and the asymptotic behavior of c a it is easy to see (at least formally) that 
for £ sufficiently large, 

Gacosier 1 ^' 3 / 4 . (B.i) 

Moreover, from standard theory that G a is smooth and positive. Biler et. al. prove in [7j that 
G a are radially non-increasing, and hence G a = G*. This same result was proved also by different 
methods in 



B.1.2 Part (ii) 

In [7], Biler et. al. show that G a is the unique, radially symmetric, self-similar solution to the PKS 
(equivalent to fixed points of the self-similar PDE (jl.lOj) and hence critical points of Q (jl.llj) ). In 
|52} [53] (see also [51]) it is shown using the moving planes method that any self-similar solution 
to (jl.ip must be radially symmetric, and hence G a is the unique self-similar solution. One can 
prove the same result with a symmetrization argument as follows. Suppose that there existed a 
non-radially symmetric self-similar solution it(£) with mass a £ (0, 8tt) and finite self-similar energy. 
Denote u(t, £) the radially symmetric solution to (ll.lOp with initial data u(0, £) = u*(£), the Riesz 
symmetric decreasing rearrangement of u. By the symmetrization inequalities in |28[ 129] . we know 
that u dominates u in the sense of mass concentration: 

u^u(t), Vt>0, (B.2) 

where if /, g are two integrable functions, f < g denotes 

/ f*(x)dx< I g*(x)dx, Vi?>0. 

J\x\<R J\x\<R 
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Since G a is the unique radially symmetric stationary point of (jl.lOp . a compactness argument 
using the energy dissipation inequality for (jl.lip implies that u{t) — > G a as r — > oo in LP for all 
p G [l,oo). Passing to the limit in (1B.2D implies 

u -< G Q . 

Further, note that any stationary solution to (jl.lOp with mass a satisfies the virial-type identity 

hence both G a and u have the same second moment. An elementary lemma regarding the Riesz 
symmetric decreasing rearrangement shows that these facts together imply u(^) = G Q (£) |34j . 

B.1.3 Part (iii) 

Follows from part (ii) using an argument similar to what is employed in [36] . 
B.1.4 Proof of (iv) 

By the above, the G a are the unique solutions to the system 

V • (G Q Vc Q ) = LG a , 
—Ac a = G a . 

The operator L can only be inverted up to the zero eigenmode, given by the Gaussian G(£), 

e H€l 2 /4 

G «> = (Sous- 

Hence, we can write G a as 

G a = qG + L _1 V • (G a Vc a ), (B.3) 

which is amenable to contraction mapping arguments. We need the following lemma regarding the 
linear operator L. 

Lemma 13. (i) L satisfies the following for all p £ (1,2] and m > 4, 

Hi-V/IU^) < ||/||lp M . (B.4) 

(ii) L satisfies the following for all p G (1,2] and m > 4, 

HVL^V/m^ll/IU+ll/llwM (B.5) 

Proof. Observe the formula, 

/•oo 

L-'Vf = - S(t- r')Vfdr'. 
Jo 

Using a compactness argument and the estimates (|2.16|) . (|2.17|) . one can justify the convergence 
of the integral for / G L g (m), q G (1,2] and in particular we see that (I2TT61) implies (fBTijh Write 
u = L _1 V/ and note that since An = V/ — • (xu) the Calderon-Zygmund inequality implies (or 
since this is L 2 , really just the Fourier transform), 

||Vn|| 2 <||/|| 2 + ||xn|| 2 <||/|| 2 + ||n|| L2(m) . 
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Hence, by C alder on- Zygmund followed by Holder's inequality and Gagliardo-Nirenberg, 



< 


\XU 4 + 


ll/IU 




< 


1 II 1 / 4 
1 "L 2 (m] 


Nlf + 11/ 


IU 


< 


1 II 1 / 4 
\ U \\L 2 (m) 




2 /2 + ll/l|4 


< 


\u\\ 1/2 


l|Vu[|^ /2 + | 


/Ik 


< 


M L 2 (m) 


+ l|V«[| 2 + 


ll/lk 


< 


\ u L 2 (m) 


+ II/II2 + II 


f\W 


< 


\ u L 2 (m) 


+ ll/lll /3 ll/ 


llf+ll/lk 


< 


\ u L 2 (m) 


+ II/II1 + II 


f||4 


< 


M L 2 (m) 


+ ll/IUp(m) 


+ II/IU, 



where the last line followed from Holder's inequality and m > 2. Hence (|B.4p implies (|B.5|) . □ 
Turning back to the proof of (iv) , we set up a contraction argument using the norm 

II/IU:=||/||l 2 M + ||v/|| 4 , (b.6) 

which by Gagliardo-Nirenberg and m > 2, embeds into every IP space for p £ [l,oo] as well as 
L q {m) for all q > 2. Indeed, recall that 

ll/Hoo<||/|l2 /3 ||V/||f <||/||x, 

and since m > 2, for p 6 [1, 2), 

ll/ll P <ll/llL 2 ( m )<H/IU- 

For p £ X, define the map Fp — > f by 

/ = aG + L- 1 V-(pVc p ), 
-Ac p = p. 

Using (IB.4j) . for any p £ (1, 2], 

ll/lli 2 (m) ^ a||G|| L 2 (m) + ||pVcp|| iP(m ). 
Using Holder's inequality, (|2.ip . and m > 2, 

||pVc p || < ||/3||L2 (m) ||Vc p ||^ 

<l|p||L 2 ( m )l|p|| P <l|p||| 2 ( m )<l|p|li- 

Hence, 

||/|| L2(m) < a||G|| L2(m) + \\pf x . 
A similar argument also shows that if fi = F[pi\, 

||/l - /2||L 2 (m) < (||Pi||l 2 (to) + l|P2||L 2 ( m )) \\Pl ~ />2||l 2 (»i)- 
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Hence for m > 2, using ()B.5p . (|2.ip and similar estimates to above, 

l|V/|| 4 < a\\VG\U + ||pV C/ ,|| 4 + ||/oVc p || LP{m) 

< Q||VG|| 4 + ||Hloo||p||4/3 + \\p\\h(m) 

< a||VG|| 4 + ||p|UI|p||i2( m ) + |Mli 2(m ) 

< a||VG|| 4 + \\p\\ 2 x , 

and similarly, 

||V/i - V/2II4 < (\\pl\\x + \\P2Wx) 1 1 Pi - P2\\x- 

Therefore, the contraction mapping theorem implies that for a sufficiently small, provided m > 2, 
^ o, which in turn proves the estimate stated in (iv). 

B.1.5 Part (v) 

The argument for (v) is similar to that for (iv). For all a E (0, 871"), G a is the unique solution of the 
system 

V • (G Q Vc a ) = LG a 
-Ac a = G a . 

Defining / := G a — Gp, subtracting the systems satisfied for each self-similar solution and rearrang- 
ing we get the following elliptic system for /: 

V • (/Vc) = (L — A a )f (B.7) 
-Ac = f 

The linear operator L — A a has a zero eigenmode with spectral projection defined by 

Pof = ( I fdx] Eq , 



where Eq is the eigenfunction in L 2 (G Q 1 c?^) corresponding to the zero mode and is given by 

E% = ^G x \ x=a . 

As above, L — A a can only be inverted up to this zero mode, and hence, 

f = P f + (L-A a )- 1 V-(fVc). 
With X defined in (IB. 6ft . one can bootstrap the result of Lemma [17] in Appendix $C]to prove, 

W\\x < a \a-P\. 

Hence, the same contraction mapping argument used to deduce (iv) can be used here provided we 
have the analogous lemma. 

Lemma 14. (i) L — A a satisfies the following for all p £ (1, 2] and m > 2, 

IKi-A^V/lk^M^II/IUpM- (B.8) 

49 



(ii) L — A a satisfies the following for all p £ (4/3, 2] and m > 4, 

||V(L - A^v/m < H/IU + ||/IUpM- (B.9) 

Proof. Inequality (|B.8P follows similarly to (|B.4|) . using the estimates collected in Proposition 
Denote u = (L — A a )~ 1 Vf as the unique solution in L^(m) to 

Lu — A a u = V/, 

and denote c the solution to —Ac = u. Similar to the proof of (|B.5P we use the Calderon-Zygmund 
inequality which implies (using also (12. IB ). 



||V«|| 2 < ll/lb + \\xu\\ 2 + \\G a Vc\\ 2 + ||t*VCa||2 
<a II/II2 + IMlL2(m) + l|Vc|| 4 + ||«|| 2 
<a ll/lb + ||w||x2(m) + IMU/3 
<a ll/lb + lk|| L 2 (m) , 

where the last line followed from m > 2. ^From here, the corresponding L 4 estimate implies (|B,9I ) 
similar to above in the proof of ()B.5p . □ 



C Appendix: Spectral Gap Estimate 

In this Appendix we sketch an independent proof of TheoremlU due to J. Campos and J. Dolbeault. 
Although this proof is independent, it was completed after knowledge of their result. In what follows 
define ^(G" 1 ^) := {/ £ I?{G~ x d£) : f fd£ = 0}. 

We prove the following result, which is equivalent to Theorem [JJ 

Proposition 13. Let f £ L^G^ 1 d^) . Then for all a G (0, 8ir) there exists some K a > such that 

HTa(r)/|| i2{G - ld0 < a e-^WfW^a^, (C.l) 

where K a and the implicit constant only depend on K < 8ir for all a < K. 

Remark 4. That Proposition [13] holds for a sufficiently small with uniformly controlled implicit 
constant and K a for a \ can be shown either by the work of |12| or by an argument essentially 
the same as Lemma W2\ using Proposition [1] (iv). Indeed, one can show K a ~ 1/2 for a small in the 
same sense as Lemma fT2l 

Let K{x) := ^log|x| be the fundamental solution for Poisson's equation in two dimensions. 
The following characterization of G a is important for what follows: 

VlogG a = Vc a -ie (C2) 
Formally linearizing Q (defined in (jl.lip ) around the stationary point G a yields (using (|(J.2|) ) . 

G(G a + ef) = g(G a ) + j 

which suggests a natural Lyapunov function for the linearized problem. Hence, define 



JZt 

G n 



d£ 



f>C * fdi 



0(e 3 
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It will turn out that F is convex and that the linear evolution is the corresponding gradient flow 
with the appropriate metric. The first step is the following dissipation inequality. The proof is a 
direct computation using (lC.2j) which we omit for brevity. 



Proposition 14. Let /o be mean zero with F(fo) < oo and let /(t) - r 

2 



G 



V ( -t- 

Gry 



Vc 



r(L-A a ) fo = Ta (r)f . Then, 
dC := -£>(/), (C.3) 



where —Ac = /. 



Drawing intuition from classical Bakry-Emery analysis and the more recent insights of entropy 
dissipation methods (see e.g. [U [201 E]) it is expected that convexity and coercivity of F are 
equivalent to the decay estimate (jC.ip . Since F is quadratic, this is in turn equivalent to strict 
positivity. 

Proposition 15. The energy F is strictly positive, that is, if f £ L^G" 1 ^) and f ^ then 
F{f) > 0. This is equivalent to strict positivity of the dissipation: 



D(J) 



G 



V ( — 

Gry 



Vc 



d* > 0. 



Proof. Linearization of Q with smooth, compactly supported perturbations and passing to the limit 
shows that necessarily F(f) > 0. Now, consider the possibility that / G Lq(G~ d£) and F(f) = 0. 
The dissipation inequality ()C.3p implies that 



Go 



V-^-Vc 



d£ = 0, 



which since G a is strictly positive implies 



V-^-Vc = 0, 

Lin 



(C.4) 



almost everywhere. A bootstrap argument using the smoothing effect of the nonlocal term shows 
that / is necessarily smooth. Moreover, we also have, 

— c = K 

for some constant K. Since / G Lq(G~ d£) we also have c G L°°, and hence /G" 1 G Taking 
the divergence of (|C.4p implies that 

Re-naming /i = /G" 1 , we see that the question of whether or not F is strictly positive reduces to 
whether or not there are any bounded, finite energy J\Vh\ 2 dx < oo, solutions to the elliptic PDE 



Ah + G a h = 0. 



(C.5) 



Ruling out bounded solutions of (|C.5|) turns out to be by far the most difficult step in the proof of 
Proposition [131 

Lemma 15. The elliptic PDE ()C,5P admits no bounded, finite energy solutions. 
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Proof. The proof requires a lengthy ODE argument and several lemmas. The first step is to consider 
any potential solution and decompose into the radial Fourier series (with £ = (r sin 0,r cos 6) T ): 

oo 

HO = £ fn{r)e ine . 

n=— oo 

Let us first rule out radially symmetric solutions f(r) := fo(r). In this case, (|C.5|) becomes 

-(rf'y + G a f = f" + -f + G a f = 0. (C.6) 
r r 

The next lemma is a standard ODE result: 

Lemma 16. The ODE (|C.6p /ias iw;o linearly independent solutions, and the possible behaviors at 
zero and infinity are f(r) ~ K for some constant K or f(r) ~ logr. 

Hence it suffices to exhibit a solution to (]C.6P which is bounded at zero and unbounded at 
infinity Such a solution will be provided by the zero eigenfunction: 

Indeed, by (ICT 

A ^ + E° a = 0. (C.7) 

By (l07l) . e(|f|) := S^G^ 1 ^) solves (EH)- Moreover, since / E%{C)dx = 1, e(r) is necessarily 
logarithmically unbounded at infinity (by (IC.5j) ), 

Lemma 17. The function e(r) is bounded at zero. 



Proof. Define, 



n(t) = / E° a (x)dx, 

J\x\<Vt 



which solves the ODE ^ 

An" + n H (nm' + nm) = 0, 

irt ' 

with boundary conditions n(0) = 0, n(oo) = 1, where m(t) := J\ x \ <y qG a (x)dx. Note that this ODE 
is linear in n{t). ^From here one can apply an analysis similar to what is done in Lemma 4.1 in [7] 
to prove that n'(0) exists and is finite, and hence E a and e are bounded at zero. The argument in 
[7] is already localized to a small neighborhood of zero, which is necessary as n(t) does not satisfy 
the same monotonicity properties that m(t) does. That m'(t) is well-behaved and satisfies certain 
monotonicity properties is necessary for the proof. Since the argument is a little technical and 
follows that of [7] very closely we omit it for brevity. □ 

Remark 5. By Lemma [T7| E^ is bounded, hence from (IC.7j) a bootstrap argument implies that 
E^ is smooth and strictly positive. 

Now we turn to the angular modes. In this case we get the ODE (which holds for the real and 
imaginary parts of the solutions) 

1 n 2 

-(rf'y--f + G a f = 0, 
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which we re-write as 

(rf'Y--f + rG a f = 0. (C.8) 
r 

Of course we have the corresponding classical ODE result: 

Lemma 18. The ODE (IC.8P has two linearly independent solutions, and the possible behaviors at 
zero and infinity are ~ r~ n or ~ r n . 

We first rule out bounded solutions supported in the mode n = 1, namely such that it is ~ r~ n 
when r goes to infinity and ~ r n when r goes to zero. For this, we use that 

A%^ + %G Q = 0, (C.9) 

and hence n\ = d r G a G~ l is a solution to (|C.8|) with n = 1. By definition ni(0) = 0. It is also 
follows from Proposition [T] that necessarily n\{r) is linearly unbounded at infinity. This rules out 
any bounded, non-zero solutions in the first angular mode. Moreover, from the monotonicity of G a: 
we get the important fact that n\{r) is strictly negative for r > 0. 

Now we confront n > 2. For this we will use second-order comparison principles against rai, 
similar to, for example, Chapter 8 in [25J. Suppose we have a bounded solution f(r) to (IC.8j) with 
n > 2. Therefore, near zero f(r) ~ r n and near infinity /(r) ~ r~ n . 

Lemma 19. Let f be a solution to (|C.8|) with n>2 which is bounded. Then f = 0. 

Proof. If / vanishes in an open neighborhood of zero then by unique continuation / = 0, therefore 
since / ~ r n near zero, we can assume / is strictly positive on some open set (replacing / by — / if 
necessary). Define x\ S (0, oo] by 

xi := sup {r > : < f(s), Vs G (0, r)} . 

If f(r) crosses zero at a finite value of r then x\ is the location of the first zero of /. If / remains 
positive for all time then x\ = oo. We will compare / to the strictly positive solution g(r) := —n\{r) 
of (|C.8j) with n = 1 on the interval (0,x\). Multiplying the ODE satisfied by g by / and vice- versa 
and then subtracting gives 

(rf')'g - (rg')'f -\{n 2 - l) gf = 0. 
Integrate now from to x\ (both sides will turn out to be integrable in the case x\ = oo): 

f X1 ( r f)'9 ~ (rg')'fdr = H - (n 2 - l) gfdr. (C.IO) 
Jo Jo r 

First notice that 

{rf'g - rg'f)' = (rf')'g - (rg')'f. 
However rf'g — rg'f is zero at r = 0, which implies (|C.10|) becomes 

/ , f Xl n 2 - 1 

lim rf'{r)g{r) — rg'{r)f(r) = / g(r) f (r)dr > 0. 

r->xi J r 

If xi < oo then f'{x\) < and f{x\) = which implies that the integral on the RHS must be equal 
to zero, which implies / = on (0,xi). By unique continuation, / = on [0, oo). If x\ = oo we 
have to first show that the RHS of (IC.lOP is integrable. This follows since g ~ r and / ~ r~ n as 
r — ¥ oo. Moreover, 

lim rf'(r)g(r) — rg'(r)f(r) = 0, 

r— >oo 

since g'(r) ~ 1 and f'(r) ~ r -3 . Hence in the case that x\ = oo, (|C.10|) still implies / = 0. □ 
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This completes the proof that (|C.5p has no non-trivial bounded solutions onR 2 . □ 

This in turn, completes the proof that F and D(f) are positive. □ 

Proposition [15] proves that F is positive, which implies F is convex since it is quadratic in /. 
The next proposition uses compactness arguments to confirm first that F is coercive and next that 
D(f) controls F. 

Proposition 16. For all a S (0, 8n) the following holds. 

(i) There is a constant C a £ (0, 1) such that for all f G L^G^d^), 

0< f fK*fdx<C a J \f\ 2 G- l dx. 

In particular, F is coercive: 

(1-C a ) J \f\ 2 G- 1 dx<2F{f), 

where 1 — C a > 0. 

(ii) There exists a constant K a > such that for all f G L^{G~ l d^), 

K a F(f) < D(f). (C.ll) 

Remark 6. By RemarkUl we do not need to worry about the behavior of these constants as a \ 0. 

First let us show why Proposition [16] completes the proof of Proposition [T3l By (jC.lip and 
(|C,3|) we get exponential decay of F, which combined with coercivity implies 

Wo ~ Hf) ~ p ^ e ~ Kat - 11/11 2 — e ' KaK 



This completes the proof of Proposition [T3l 

Now let us turn to the proof of Proposition [TBI 



Proposition [TR Since F(f) > we already have 

JfK*fd£< J l/l'G- 1 ^, 
for all / 7^ 0. Suppose there exists a sequence of {fkYkLii fk 7^ 0, such that 

lim / f k fC * f k dH = lim / l/^C" 1 ^. 
By homogeneity, we can assume without loss of generality that J \ f k \ 2 G~^d£, = 1- Hence, 



lim / / fc /C * f k d£ = 1. 



By the boundedness of f k in L^{G a l d^) we may extract a subsequence (not relabeled) that weakly 
converges in Lq(G~ d£) to some limit h which by lower semicontinuity satisfies 
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Since f k can also be chosen to converge weakly in L , we have that J hd£, = 0. We claim that 

1 = lim / f k JC * f k d£ = I hJC* hd£, (C.12) 

fc->oo J J 

which implies both that h ^ and F{h) < 0, in contradiction with Proposition [15] We now prove 
(1CTT21) . First, define Kf = K{\x - y\) 

^-5<\x-y\<R anci break up the convolution into 
/ fk>C*fkdx= I fkKf-*f k dx + ^- j f k {x) f k {y) \og\x -y\dxdy 

J J 47T J\ x -y\ <S 

+ ~T I fk(x)fk(y) log \x - y\ dxdy 

47r J\x-y\>R 

= Tl + T2 + T3. 

The f k are all uniformly well localized, hence the term T3 can be made arbitrarily small by choosing 
R > 1 large: 

f fk{x)f k (y) log \x-y\ dxdy < -^|^ / \fk(x)\\fk(y)\\x - y\dxdy 

J\x-y\>R K J\x-y\>R 

<^ / \fk(x)\\fk(y)\(\x\ + \y\)dxdy 

K J\x-y\>R 

S^IIAIM/MM Ik. 



<^/lAI 2 e;^. 



Similarly, since the f k are uniformly bounded in L 2 , the first term T2 can be made arbitrarily small 
by choosing 5 small: 

/ fk(x)fk(y) log \x - y\ dxdy < \\f k \\ 2 L 2 || log \x - y\ l| a ._ 1/ |< ( 5||x,i. 

J\x-y\<8 

Notice that the exact same arguments apply to h. Hence, it suffices to prove that for all 5, R we 
have 



J fkK? * f k dx -> J hKf * hdx. 



This convergence follows from classical weak convergence arguments, but we will still give a proof. 
First, note that this would be straightforward if f k converged to h strongly in L^G" 1 ^), but weak 
convergence suffices. Indeed, consider 

J fkK* * fkdx - J hKf *hdx = J (f k - h)Kf *hdx + J f k Kf * (f k - h). 

Since — /i in L 1 the first term converges to zero, so we need only focus on the latter. Define 

v k (x) := / (f k (y) ~ Ky))Kf{x - y)dy = (f k - h) * K$. 



By weak convergence again, v k (x) — > pointwise a.e.. Let v(x) := sup k v k (x). Firstly, 

\v{x)\ < sup - h\\ L i(B R ( x ^max(\log6\ ,|logi?|). 

k 



55 



Since fk and h are exponentially localized in L 1 uniformly in k, it follows that v{x) is pointwise 
exponentially localized as well as bounded, and hence integrable. By the dominated convergence 
theorem, it follows that Vk — > strongly in L l . By interpolation this implies v k — > in LP for all 
p € [1, oo) from which it follows that J fk^f * (fk — h) — > 0. Putting all of the estimates together, 
we deduce that 

J fkK- * fkdx — > J hK, * hdx, 

which is the desired contradiction. This completes the proof of (i). 

The proof of (ii) continues in a similar fashion. We suppose there exists a sequence {/fcl^Lx of 
functions fk ^ normalized such that F(fk) = 1 but D(fk) —> and derive a contradiction. By the 
coercivity estimate (i), this implies 2 < ||/^|| 2 2 x < 2(1 — C a )~ l and ||V/C*/ fc || 2 < 1. Extracting 

a subsequence if necessary, we can assume that fk converges weakly to some h in L 2 (G~ 1 d^). Since 
fk can also be chosen to converge weakly in L 1 , we have that J hd£ = and hence h 6 L^G^ 1 d£) . 
By lower semicontinuity of D(f) with respect to weak convergence, we also deduce that D(h) = 0. 
To get a contradiction, we have just to prove that h ^ 0. 



It follows from the bound ||V/C * fk\\2 i$ 1 and the boundedness of D(fk) that j G c 



2 



is uniformly bounded. Define V = V(G a ) = {F : \\F\§ = j G Q (|VF| 2 + \F\ 2 )d£, < oo}. Hence, 
V is compactly embedded in L 2 {G a d£) and since -ff- is bounded in V, it is relatively compact 
in L 2 (G a d^). Hence, fk is relatively compact in L 2 (G~ 1 (i^) and extracting a subsequence, we 
deduce that fk converges strongly to h in L 2 (G~ l d^) and hence h ^ 0. This ends the proof of the 
Proposition. 

□ 
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